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ABSTRACT 


A  three  -  component  Laser  Doppler  Velocimetry  (LDV)  system  for 
the  David  Taylor  Naval  Ship  Research  and  Development  Center's  towing 
carriages  became  operational  in  mid  1984.  This  is  one  in  a  series  of 
reports  to  help  experimenters  efficiently  use  the  new  instrumenta¬ 
tion.  The  series  also  serves  as  a  repository  for  a  large  volume  of 
information  gathered  while  purchasing,  testing,  and  optimizing  this 
velocity  measurement  system.  This  particular  report  deals  with  mean 
velocity  measurement  errors:  both  minimizing  them  and  quantifying 
what  remains. 

Thirteen  measurement  errors  are  quantified  in  terms  of: 
velocity  field  characteristics,  LDV  equipment  design,  and  user 
chosen  system  operating  parameters.  Operating  parameter  effects  on 
measurement  error  are  emphasized  because  usually  they  are  the  only 
way  an  experimenter  can  control  errors. 

Background  sections  and  appendixes  explain  in  considerable 
detail,  the  origin  of  these  measurement  errors.  First,  the  operating 
parameter  choices  for  the  towing  tank  LDV  system  are  listed  with  a 
description  of  how  they  affect  measurements.  Then,  each  of  the  meas¬ 
urement  errors  is  discussed,  and  this  helps  establish  a  functional 
relationship  to  operating  parameters,  LDV  equipment  design,  and  flow 
field  characteristics.  Various  appendixes  support  and  amplify  the 
information  given. 

Three  user  reference  sections  are  included  in  this  report.  The 
first  provides  guidance  in  making  operating  parameter  choices  and 
measurement  -  error  estimations  prior  to  an  experiment.  The  second 
describes  how,  during  an  experiment,  the  user  can  determine  if  proper 
parameter  choices  have  been  made.  The  last  section  details  how  final 
error  estimates  can  be  made  after  an  experiment,  for  inclusion  in 
written  reports. 
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ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  done  at  the  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  under  a 
project  for  Laser  Doppler  Velocimeter  System  Development.  The 
project  was  funded  by  the  6.2  Ship  and  Submarine  Technology  Program, 
Program  Element  62543N,  Task  Area  SF  43421,  and  DTNSRDC  Work  Unit 
1506-253  (FY-84)  and  1506-130  (FY-85). 

INTRODUCTION 

In  the  last  quarter  of  fiscal  year  1981  (FY-81),  work  began  on  a 
three-component  Laser  Doppler  Velocimeter  (LDV)  system  for  all  towing 
basins  at  DTNSRDC.  It  was  envisioned  that  the  primary  use  of  the 
system  would  be  to  do  ship  wake  surveys  and  measure  mean  velocities. 
A  contract  was  awarded*  to  TSI  Inc.  to  build  the  optics  system. 

Appendix  1  details  the  optical  configuration.  The  optics  system 
was  completed  and  delivered  to  DTNSRDC  in  June  1982  for  inspection 
and  checkout.  The  system  met  all  contract  specifications. 

A  new  contract  was  awarded**  to  TSI  Inc.  to  develop  a  traversing 
system  and  data  analysis  system  for  the  Towing  Tank  LDV. 

The  completed  system  was  delivered  and  installed  in  July  1983. 
However,  several  performance  specifications  were  not  met.  TSI  took 
the  system  back,  made  modifications,  and  redelivered  the  system  in 
March  1984  (Figures  1  and  2).  This  time  the  system  was  accepted. 

The  procedures  and  skills  of  aligning  and  setting  up  the  system 
are  not  covered  here.  This  information  is  available  from: 

1.  Previous  system  users 

■'  *  *  * 

2.  DTNSRDC  Departmental  Report^ 

3.  TSI  Inc.  manuals^"^ 


Contract  number  N00167-81-C-0284 ,  Oct.,  1981. 

Contract  number  N00167-82-C-0212 ,  Oct.,  1982. 

A  complete  listing  of  references  is  given  on  page  94 


There  is  essentially  only  one  way  to  physically  set  up  the  system. 
However,  once  the  LDV  system  is  in  place,  users  are  confronted  with 
many  operating  parameter  choices  before  taking  any  data. 

This  report  will  guide  users  in  making  operating  parameter 
choices  as  well  as  enable  them  to  quantify  the  errors  present  in 
their  measurements.  Capabilities  and  errors  are  a  function  of  the 
LDV  equipment  design  and  flow  field  characteristics.  However  they 
are  also  highly  dependent  on  operating  parameters  chosen  by  users. 

The  sections  of  this  report  are  comprised  of  two  types: 

1.  Background  Information 

2.  Reference  Information 

The  following  two  sections  supply  background  information.  The  first 
lists  the  operating  parameter  choices  for  the  towing  tank  system  and 
describes  how  they  affect  LDV  measurements.  The  second  lists  possi¬ 
ble  measurement  errors  and  establishes  their  functional  relationship 
to  operating  parameters,  LDV  equipment  design,  and  flow  field 
characteristics.  Various  appendixes  support  and  amplify  the  informa¬ 
tion  given. 

Three  reference  sections  follow.  The  first  provides  guidance  in 
making  operating  parameter  choices  and  measurement  error  estimations 
pr i or  to  an  experiment.  The  second  describes  how,  during  an  experi¬ 
ment,  the  user  can  determine  if  proper  parameter  choices  have  been 
made.  The  last  section  details  how  final  error  estimates  can  be  made 


LDV  OPERATING  PARAMETERS 

PRINCIPLES  OF  LDV  MEASUREMENTS 

Before  the  choice  of  various  operating  parameters  can  be  made, 

the  basics  of  the  LDV  measurement  technique  must  be  understood.  The 

principles  behind  dual  beam,  LDV  measurements  are  described  in  great 

p  Q  1  n 

detail  in  many  publications  The  technique  description  that 

follows  will  be  much  briefer  and  less  detailed. 

Measurement  Basics 

LDV  systems  typically  focus  and  cross  a  number  of  laser  beams  at 
a  measurement  "point".  The  towing  tank  system  has  five  separate 
beams  of  two  colors.  Particles  in  the  water,  passing  through  the 
beam  crossing  region  scatter  light  that  the  LDV  system  collects  and 
utilizes  in  determining  the  particle  and  fluid  velocity  (if  the 
particle  "follows"  the  flow). 

The  determination  of  particle  velocity  is  best  understood  in 
terms  of  interference  "fringes"  of  light  effectively  set  up  within 
the  beam  crossing  region.  The  fringes  are  a  series  of  parallel, 
alternating  light  and  dark  planes.  The  intensity  of  the  scattered 
light  changes  as  the  particle  moves  between  light  and  dark  fringes. 
The  fringes  have  a  uniform  and  accurately  determinable  spacing 
(fringe  spacing,  S,  is  a  function  of  laser  beam  wavelength  and  beam 
crossing  angles).  A  photomultiplier  and  an  electronic  signal  proces¬ 
sor  determine  the  time  required  by  the  particle  to  cross  a  given 
number  of  fringes  (Figure  3).  This  time  is  passed  to  the  LDV 
.'ysfem's  computer  for  storage.  Later,  the  computer  is  used  to 
calculate  a  component  of  the  particle’s  velocity  by  multipling  the 
known  fringe  spacing  and  the  fringe  crossing  frequency  (obtained  from 
stot.  a  rime  data).  By  looking  at  scattered  light  from  different  sets 
of  beams,  all  three  velocity  components  are  determined. 

Simple,  stationary  fringes  allow  determination  of  velocity 
component  magnitudes.  But  there  is  still  a  need  to  distinguish 
b- ■  f  wro  particles  traveling  at  the  same  speed  but  in  opposite 


directions.  If  fringes  could  move  and  move  fast  enough,  they  would 
always  sweep  across  particles  in  the  same  direction,  regardless  of 
the  particle's  direction  of  motion.  A  fringe  speed  greater  than  the 
fastest  particle  moving  in  the  same  direction  would  be  required. 

The  towing  tank  LDV  system  effectively  allows  fringes  for  each 
component  to  move  past  a  fixed  point  at  frequencies  (called  a 
frequency  "shift"  or  FS)  between  .002  and  10  MHz,  inclusive.  This 
frequency  is  precisely  the  frequency  of  scattered  light  intensity 
changes  that  the  LDV  system  will  measure  for  a  stationary  particle. 
Particles  moving  against  the  fringes  will  have  higher  frequencies. 
Particles  moving  in  the  same  direction  as  the  fringes  will  have  lower 
frequencies.  In  data  processing,  the  shift  frequency  must  be 
subtracted  from  the  measured  frequency  before  multiplication  by  the 
fringe  spacing  yields  the  particle  velocity.  Any  resulting  negative 
velocities  indicate  particles  moving  with  the  fringes,  while  positive 
velocities  indicate  motion  in  the  opposite  direction. 

Spatial  Resolution 

The  LDV  measurement  signals  only  occur  when  a  particle 
simultaneously  scatters  light  from  at  least  two  laser  beams.  Thus, 
it  is  apparent  that  LDV  measurements  are  particle  velocities  that  oc¬ 
cur  somewhere  within  the  laser-beam  crossing  region  (henceforth 
called  the  "measurement  volume") .  Spatial  resolution  is  at  least  as 
good  as  these  measurement  volume  dimensions*.  It  should  be  noted 
that  the  scattered  light  from  different  sets  of  beams  are  used  in 
determining  each  of  the  three  components  of  velocity.  Figure  4 
indicates  which  beams  are  involved  and  the  measurement  volume  dimen¬ 
sions  for  each  component. 

Laser  beam  cross-sections  have  no  sharply  defined  edges.  One 


Some  LDV  systems  have  scattered  light  receiving  optics 
focused  only  on  portions  of  the  beam  crossing  region  and  there¬ 
fore  get  better  spatial  resolution. 


II 


beam  diameter  definition,  "Dg_2»”  is  based  on  the  circle  of  points 
where  the  light  is  13.5%  of  the  intensity  found  at  the  beam  center. 
Beams  with  a  circular  cross-section  of  diameter  D^_2  are  used  in  the 
geometric  calculation  of  the  measurement  volume  dimensions  (Figure 
4) .  A  more  meaningful  definition  of  measurement  volume  boundaries  is 
where  some  average  sized  particle  has  a  scattered  light  signal  just 
above  detection  thresholds. 


Measurement  Response  Time 

The  signal  processors  and  computer  of  the  towing  tank  LDV  system 
determine  the  particle  and  water  velocity  on  the  basis  of  timing  "M^" 
fringe  crossings.  Velocity  changes  that  occur  during  the  crossing  of 
fringes  cannot  be  detected  but  are  averaged  together. 

"^avg  =  Mf*Sy/(FS*Su+  U)  (1) 

t  1  ♦  *  * 

where  U;  measured  velocity  component 
Sy ;  component  fringe  spacing 

Note  that  this  averaging  time  is  a  function  of  both  particle  (or 
water)  velocity  and  the  operating  parameter  of  frequency  shift. 

The  signal  processors  also  have  an  "electronic  downtime"  after 
each  measurement  due  to  the  fringe  timing  circuits  and  the  transmis¬ 
sion  of  data  to  the  computer.  The  manufacturer,  TSI  Inc.  lists  this 
downtime  as; 


^dl  =  8  MS  (2) 

velocity  measurements  cannot  be  made  at  time  intervals  less 

T 
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response  of  the  LDV  system  is  limited  also  by  the  type  of 


Sepa  rate 
than 

The 


I 

i 

f: 


iv"! 


i 

f:-' 

L  • 


i 


i^. 


i'- 

1 

K' 

C'' 


^i^i^w^^yy^y^ij^MBfiu»nnmjin.im.ii»»)i.ii 


signal  processors  used  in  the  towing  tank  system.  These  counter  type 
processors  work  best  when  only  one  particle  is  ia  the  measurement 
volume  at  a  time.  Randomly  arriving  particles  can  only  satisfy  this 
condition  if  the  measurement  volume  is  frequently  void  of  particles 
for  periods  longer  than: 

’'d2  =  Omv  /  W^+v2)»-5  ,3) 

where  ^mv'  short  axis  dimension  of  LDV  Meas.  Vol. 

U,V:  velocity  components  (see  Figure  4) 

This  period  might  be  called  a  particle  "d'^wntime." 

Velocity  fluctuation  frequencies  greater  than  1/(2*T^2)  are  not 
detectable  by  the  LDV  system  under  even  optimal  particle  arrival 
rates.  The  towing  tank  system  usually  sees  much  slower  particle  ar¬ 
rival  rates  that  greatly  reduce  the  upper  frequency  limit  (50  Hz  or 
less)  on  detectable  velocity  fluctuations.  However,  if  the  fluctua¬ 
tions  are  repetitive  (such  as  those  due  to  a  propeller  rotating  at 
constant  RPM)  velocity,  change  frequencies  of  1/  {2*’V  are  obtain¬ 

able  for  data  collected  and  averaged  over  repeated  fluctuation  cy¬ 
cles. 

FREQUENCY  SHIFT  OPERATING  PARAMETERS 

This  section  describes  the  operating  parameter  choices  allowed 
by  the  frequency  shift  components  of  the  towing  tank  LDV  system. 
These  parameters  are  set  by  the  orientation  of  Bragg  cells  in  the  LDV 
optico  system  and  by  push  button  controls  on  the  accompanying 
electronic  control  boxes.  The  design  and  operation  of  the  Bragg 
cells  and  their  signal  mixing  electronics  is  covered  extensively  in 
TSI  manuals.^'^ 

Different  frequency  shifts  are  not  obtained  by  actually  changing 
the  speed  of  the  moving  fringes.  The  LDV  optics  system  can  use  Bragg 
cells  operating  at  only  a  single,  unadjustable  frequency.  However, 
the  LDV  photomultiplier  output  signal  (fixed  frequency  shift)  can  be 
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electronically  mixed  with  a  signal  of  different  and  adjustable 
frequency.  The  voltage  signal  obtained  behaves  like  an  LDV  signal 
with  a  frequency  shift  equal  to  the  difference  between  the  Bragg  cell 
and  mixing  signal  frequencies  This  derived  LDV  signal,  with  its 
"effective"  frequency  shift,  is  used  by  the  towing  tank  LDV  signal 
processors . 

Effective  frequency  shifts  are  available  in  the  range  of  .002  to 
10  MHz.  Obtainable  values  are  in  discrete  steps  of 
. 002 , . 00 5 , . 010 , . 02 , . 05 , . 10 , . . . .  The  effective  direction  of  fringe 
motion  is  also  selectable.  One  additional  fringe  motion  frequency  is 
available  for  each  velocity  component  when  the  mixing  circuitry  is 
bypassed  (100  MHz  for  the  streamwise  component,  40  MHz  for  the  verti¬ 
cal  component,  20  MHz  for  the  on-axis  component).  The  direction  is 
unchangeable  in  this  case. 

The  frequency  shift  for  each  velocity  component  is  chosen  ac¬ 
cording  to  the  range  of  expected  velocities.  The  lowest  frequency 
shift  that  allows  M^  fringe  crossings  for  these  veloci’.if.s  i the 
choice  that  will  usually  minimize  overall  measurement  error. 

COUNTER  PROCESSOR  OPERATING  PARAMETERS 

The  counter  processors  for  the  towing  tank  LDV  are  manufactured 
by  TSI  Inc.  (Model  1990).  There  is  a  separate  processor  for  each  of 
the  three  velocity  component  channels.  The  LDV  signal  input  into  the 
processor  is  a  voltage  that  varies  almost  sinusoidally  (under 
noise-free  conditions).  Positive  voltage  peaks  occur  when  the  parti¬ 
cle  is  passing  through  a  light  fringe  and  negative  voltage  peaks  oc¬ 
cur  when  a  dark  fringe  is  crossed.  The  magnitudes  of  the  peaks  are 
modulated  by  the  position  of  the  particle  in  the  measurement  volume. 
The  largest  magnitudes  occur  when  the  particle  is  illuminated  most 
intensely,  near  the  center  of  the  measurement  volume.  A  typical  LDV 
signal  (called  a  signal  "burst")  input  to  the  processors  is  shown  in 
Figure  3.  The  processor  output  is  a  digital  signal  containing  the 
count  of  clock  cycles  (1  ns  each)  occurring  while  the  particle 
crosses  M^  fringes. 
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The  counter  processors  have  many  different  parameters  that  must 
be  set.  They  can  basically  be  divided  into  those  that  change  the 
input  signal  and  those  that  affect  the  fringe  counting  method.  A 
very  complete  description  of  the  processor  and  parameter  settings  is 
found  in  Reference  3.  This  section  will  provide  only  the  details 
needed  for  later  estimation  of  velocity  measurement  errors. 

Signal  Conditioning 

A  selected  high  pass  and  low  pass  filter  are  first  applied  to 
the  LDV  signal.  Twelve  settings  are  possible  for  each  of  the  high 
and  low  pass  filters^.  Each  setting  has  a  minimum  roll  off  of  30 
dB/octave.  These  filters  should  be  set  to  eliminate  as  much  signal 
noise  as  possible  and,  yet,  not  filter  out  any  of  the  expected  signal 
frequencies  . 

A  continuously  variable  gain  control  (-31  dB  to  +34  dB)  is 
available.  Signal  voltage  amplitudes  must  exceed  counter  thresholds 
of  +/-  50  mV  to  trigger  the  timing  circuitry.  The  upper  limit  on  an 
acceptable  gain  setting  is  controlled  by  the  necessity  of  keeping 
signal  noise  amplitudes  well  below  the  counter  thresholds  (to  avoid 
spurious  counts).  Smaller  signal  gains  and,  hence,  noise  amplitudes 
reduce  errors  in  the  processor's  timing  of  fringe  crossings. 

Larger  particles  tend  to  produce  more  scattered  light  and, 
therefore,  higher  voltage  input  signals.  An  amplitude  limit  control 
throws  away  processor  timings  done  on  signals  exceeding  a  maximum 
voltage.  A  good  choice  of  this  maximum  voltage  parameter  can 
eliminate  data  from  large  particles  not  following  the  water  flow. 

Processor  Operating  Method 

Three  types  of  timings  can  be  done  by  the  processor  on  the  LDV 
sionil.  A  "primary  fringe  timing"  counts  the  processor  clock  cycles 
ovc’r  a  set  number,  M^,  of  fringe  crossings.  A  "secondary  fringe  tim¬ 
ing"  counts  the  clock  cycles  over  a  smaller  number  of  fringe  cross¬ 
ings.  The  primary  and  secondary  fringe  timings  occur  simultaneously. 
Finally,  a  "burst  timing"  counts  clock  cycles  over  all  of  the  fringe 


crossings  that  a  single  particle  makes  as  it  crosses  the  measurement 

volume  . 

There  are  four  "modes"  of  processor  operation  that  can  be 
chosen.  Each  uses  the  three  timings  in  a  different  way.  In  continu¬ 
ous  (Cont>  mode,  as  many  sets  of  primary  and  secondary  fringe  timings 
as  possible  are  made  on  each  particle  (or  signal  burst).  In  each 
set,  primary  and  secondary  timings  are  compared  for  consistency.  If 
consistent,  the  primary  fringe  timing  and  fringes  counted  (M^)  are 
the  outputs. 

In  the  single  measurement  per  burst  (SM/B)  mode,  only  one 
primary  and  secondary  fringe  timing  is  made  for  each  particle.  They 
ate  compared  for  consistency.  The  single  primary  fringe  timing  and 
the  fringes  counted  (M^)  are  the  outputs. 

Total  burst  count  (TBC)  mode  also  has  an  output  of  a  single 
primary  fringe  timing  (if  consistent  with  the  secondary  fringe 
timing);  but  instead  of  the  total  number  of  fringes  encountered 

(as  used  in  burst  timing)  is  the  output. 

Finally,  in  the  total  burst  mode  (TBM)  the  burst  timing  and  the 
total  number  of  burst  fringe  crossings  (if  less  than  256)  are  the 
outputs.  There  is  no  consistency  check. 

Separate  from  the  mode  selection  switch  is  a  fringe  count  set¬ 
ting.  This  sets  the  value  for  primary  fringe  timing  in  the  Cont, 
SM/B,  and  TBC  inodes.  Values  of  2,  4,  8,  16,  32,  64  ,  or  128  can  be 
selected.  Secondary  timings  are  done  for  1/2  fringe  crossings 

when  ig  less  than  8.  Otherwise  5/8  fringe  crossings  are  used. 

In  the  TBM  mode  the  fringe  count  setting  is  the  minimum  number  of 
fringe  crossings  for  the  burst  timing  measurement  to  be  accepted  as 
v,i  lid. 

1  "r-omparison  %"  control  sets  the  allowable  velocity  difference 
vn‘>n  t.r  imary  ind  secondary  timings  are  checked  for  consistency. 
Grrnter  differences  cause  the  data  to  be  thrown  out.  A  comparison  % 
'  1  ue  of  0  I  indicating  no  consistency  check)  must  be  chosen  for  the 
TBM  node  where  no  consistency  check  is  possible. 


maximum  timing  accuracy.  Processor  timings  consist  of  a  clock  cycle 
count  digitally  stored  as  a  12-bit  mantissa  and  a  4-bit  exponent. 
The  automatic  setting  insures  that  the  mantissa  count  is  as  high  as 
possible  (and  the  exponent  count  is  as  low  as  possible) .  Conceivable 
timing  errors  of  +/-  1  in  the  mantissa  then  have  the  least  possible 
effect  on  measurement  accuracy. 


COMPUTER  SYSTEM  OPERATING  PARAMETERS 

This  section  describes  the  operating  parameter  choices 
encountered  when  transforming  counter  processor  outputs  into  esti¬ 
mates  of  the  three  components  of  mean  velocity.  Some  parameters  deal 
with  which  and  how  many  of  the  three  counter  processor  outputs  are 
accepted  and  stored  (i.e.,  data  acquisition).  Other  choices  concern 
whether  data  corrections  are  made  based  on  experimental  calibrations 
(i.e.,  data  correction).  The  rest  of  the  parameters  in  this  section 
deal  with  how  the  stored  data  is  later  manipulated  to  give  the  best 
possible  estimate  of  the  mean  velocity  (i.e.,  data  manipulation). 

Reference  4  details  parameter  selection  and  the  operation  of  the 
DEC  computer  interface  used  in  the  towing  tank  LDV.  Reference  7 
gives  a  brief,  global  description  of  standard  TSI  software  ("DRP-3") 
for  data  acquisition  and  manipulation.  The  software  for  the  towing 
tank  system  is  a  modified  version  of  DRP-3. 


Data  Acquisition 

The  control  of  the  data  flow  from  the  counter  processors  to 
computer  storage  is  done  by  Interface  Modules  (TSI  Inc.,  Model  1998) 
and  a  direct  memory  access  (DMA)  board.  This  interface  can  operate 
in  the  random  or  in  the  coincidence  mode.  In  random  mode,  data  from 
any  processor  is  sent  as  soon  as  it  is  ready.  All  processible  sig¬ 
nals  received  by  each  counter  are  stored  except  those  that  are  ready 
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during  the  time  of  data  transmission  from  a  processor  (a  period  of  8 
* 

r\s  )  .  Coincidence  mode  data  storage  requires  valid  data  points  from 
all  three  counter  processors  at  the  "same  time".  The  "same  time"  is 
arrival  within  a  coincidence  time  window  that  can  be  set  anywhere 
between  10  qs  and  10  ms.  The  time  window  starts  with  the  first  valid 
data  point  from  any  processor.  If  data  from  the  other  two  processors 

arrives  within  the  time  window  period,  all  three  data  points  are 

,  _  ★  ★ 

transmitted  and  stored. 

A  computer  interface  option  allows  transfer  of  the  time  between 
measurements.  An  additional  16-bit  word  is  sent  with  each  measure¬ 
ment.  It  contains  (when  properly  decoded)  the  time  delay  between 
valid  measurements. 

The  final  operating  parameter  in  data  acquisition  is  the  number 
of  valid  measurements  accepted  at  each  measurement  location.  The  sum 
of  the  contributions  from  each  processor  (or  velocity  component  chan¬ 
nel)  can  be  set  in  computer  software.  Operation  in  a  purely  random 
transfer  mode  usually  yields  more  measurements  of  some  components 
than  others.  However,  a  device  in  the  interface  insures  that  256, 
512,  1024,  or  2048  data  points  from  each  channel  are  transferred.  An 
equal  number  of  measurements  in  each  channel  (i.e.,  for  each  velocity 
component)  is  also  assured  if  operation  is  in  the  coincidence  mode. 
However,  in  this  mode,  finite  -  time  windows  mean  some  data  points 
will  be  thrown  away  because  of  noncoincidence.  Therefore,  the  rate 
of  data  transfer  to  the  computer  will  be  lower. 


Data  Correction 

The  collected  data  should  usually  be  corrected  in  one  or  more 
ways  based  on  calibration  tests  performed  after  the  LDV  system  is  in 
place.  A  few  measurements  taken  before  a  test  begins  (and  after 


*  Reflects  instantaneous  storage  rate  during  the  filling  of 
certain  finite  buffers  (32K  16-bit  words) 

**  A  required  processor  order  can  also  be  set.'^ 
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every  subsequent  LDV  optics  realignment)  can  yield  corrections  for 
the  effects  of; 


1.  Velocity  component  directional  alignment 

2.  LDV  probe  strut  flow  disturbance 

3.  Frequency  shift  and  signal  mixing  accuracy 

The  LDV  user  can  choose  to  later  adjust  his  experimental  results  us¬ 
ing  these  corrections. 

Experience  has  shown  that  checking  the  directional  alignment  of 
the  on-axis  component  is  especially  important.  Even  very  careful 
optical  alignment  may  allow  directional  errors  of  as  much  as  3°. 

It  is  desirable  to  correct  mean  velocity  estimates  for  the  flow 
disturbance  of  the  LDV  probe  strut  on  the  ship  model  flow  field  (Fig¬ 
ure  2) .  Undisturbed  velocity  values  are  always  sought  in  towing  tank 
tests.  The  actual  flow  disturbance  will  be  different  for  every 
experiment  and  for  every  measurement  location.  However,  a  first 
order  correction  is  possible.  It  is  based  on  the  measured  velocity 
disturbance  induced  at  the  measurement  volume  by  towing  only  the  LDV 
probe  strut. 

The  effective  frequency  shift  may  not  always  be  exactly  as 
indicated  by  the  shift  controller.  However  in  calculating 
velocities,  the  computer  software  uses  the  indicated  shift  exactly. 
This  results  in  possible  errors  if  no  correction  is  made. 

All  of  the  above  corrections  are  possible  from  calibration 
velocities  taken  when  the  probe  strut  is  towed  alone.  Then  any 
non-zero,  mean  velocity  (in  the  vertical  or  on-axis  direction)  is  due 
to  these  uncorrected  effects.  Streamwise  mean  velocity  differences 
from  the  tow  speed  are  also  due  to  these  same  effects  along  with 
fringe  spacing  measurement  error. 

Data  Manipulation 

The  straight-forward  arithmetic  mean  of  measurements  is  not 
always  the  best  way  to  get  accurate  mean  velocity  estimates  from 
counter  processor  data.  Because  processor  measurements  are  not 
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random  samples  of  the  water  velocity,  certain  data  biases  affect 
arithmetic  means  of  the  data. 

Processor  measurements  occur  only  when  "measureable"  (a  minimum 
number  of  strong  fringe  crossings  must  be  timed  for  a  valid  measure¬ 
ment)  particles  arrive  at  the  measurement  volume.  The  rate  of  parti¬ 
cle  arrival  is  a  function  of  the  velocity  direction  and  magnitude 
(velocity  bias).  The  "measureabil i ty"  of  any  particle  is  a  function 
of  the  particle  velocity  direction  (fringe  bias).  Selected  processor 
frequency  filters  may  attenuate  and  reduce  "measureabil ity"  of  some 
valid  LDV  signal  frequencies  more  than  others  (filter  bias). 

The  LDV  system  computer  software  contains  an  option  that  through 
data  manipulation  can  correct  for  velocity  bias.  The  key  to  the  cor¬ 
rection  is  information  on  the  measurement  volume  transit  time  for 
each  particle.  Thus,  this  correction  can  be  done  only  if  the  proces¬ 
sors  are  operated  in  TBC  or  TBM  modes.  Another  method  of  bias 
correction^^  can  eliminate  both  velocity  and  fringe  bias  if  all  three 
velocity  components  are  measured  coincidentally  ( i  .e .  within  a  time 
window  shorter  than  that  needed  for  significant  velocity  changes) . 
Also,  both  fringe  and  filter  bias  can  effectively  be  eliminated  by 
proper  choice  of  operating  parameters  already  mentioned  (particularly 
the  channel  frequency  shift  and  counter  processor  input  filter  set¬ 
tings)  . 

Finally,  when  large  amounts  of  LDV  measurements  are  taken,  a 
small  number  of  invalid  measurements  are  sometimes  also  stored.  Usu¬ 
ally  these  invalid  measurements  (often  caused  by  signal  noise)  are 
detectable  because  of  their  marked  disparity  with  the  overwhelming 
majority  of  the  rest  of  the  measurements.  Measurement  error  analysis 

often  assumes  a  Gaussian  distribution  of  measurements  about  a  mean 

1  2 

and  then  applies  "Chauvenets'  Criterion"  .  This  eliminates  data 
values  that  are  too  many  standard  deviations  from  the  data  mean  and 
therefore  unlikely  to  be  part  of  the  distribution  of  data.  The 
computer  software  has  data  manipulation  routines  that  allow  such 
"invalid"  data  values  to  be  found  and  eliminated. 
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OPERATING  PARAMETERS  SUMMARY 


The  operating  parameter  choices  described  in  this  section  have 
come  under  three  headings.  Table  1  lists  all  the  parameters  under 
their  appropriate  heading  for  easy  refeiance.  Certain  parameter 
choices  should  always  be  the  same  for  towing  tank  LDV  system  setups. 
These  choices  are  listed  in  Table  1  and  will  be  assumed  in  the  sec¬ 
tions  to  follow. 


TABLE  1  -  TOWING  TANK  LDV  SYSTEM  PARAMETER  CHOICES 


FREQUENCY  SHIFT  OPERATING  PARAMETERS: 

I- A  Frequency  Shift . 

COUNTER  PROCESSOR  OPERATING  PARAMETERS: 

II- A  Input  Frequency  Filters . 

II-B  Input  Gain . 

II-C  Amplitude  Limit . 

II-D  Mode  of  Operation . 

II-E  Mf  or  Fringe  Crossings  Counted 

II-F  Count  Comparison  Accuracy . 


I I - G  Timer  Exponent  Control  . Automatic 

COMPUTER  SYSTEM  OPERATING  PARAMETERS: 

III- A  Random  or  Coincident  Transmission .  Random 

III-B  Time  Between  Measurements . Not  Stored 

1 1 1 -C  Total  Measurements  per  Location . 

III-D  Distribution  of  Measurement  Among  Components .  Equal 


III-E  Velocity  Component  Direction  Correction..,, 
III-F  LDV  Probe  Strut  Flow  Disturbance  Correction 
1 1 1 -G  Frequency  Shift  -  Signal  Mixer  Correction.. 
1 1 1 -H  Velocity  Bias  Correction.... . 


III-I  Fringe  Bias  Correction . . .  None 

III-J  Filter  Bias  Correction .  None 


III-K  Erroneous  Data  Elimination 


VELOCITY  MEASUREMENT  ERRORS 

In  this  section,  contributors  to  towing  tank  LDV  velocity  meas¬ 
urement  error  will  be  examined  and  quantified.  They  are  divided  into 
three  groups.  The  first  group  of  errors  affect  the  accuracy  of 
individual  particle  velocity  component  measurements.  The  second 
group  of  errors  affect  the  accuracy  of  transforming  those  particle 
measurements  into  fluid  mean  velocity  estimates.  The  third  group 
enters  through  the  manipulations  necessary  to  obtain  final  (with  cor¬ 
rections)  ,  dimensionless  (based  on  towing  speed),  undisturbed  (by  LDV 
probe  strut),  mean  velocities. 

Error  magnitudes  cannot  be  discussed  as  if  they  have  a  certain 
value.  Rather  sets  of  individual  measurements  and  sets  of  mean  meas¬ 
urements  will  have  a  distribution  of  values  about  the  true  mean 
value.  Error  magnitudes  ([measurement  -  true  mean  value |)  have  a 
related  distrubition.  This  section  will  characterize  velocity- 
component  error  distributions  by  calculating  their  standard  devia¬ 
tions  (nondimensionalized  by  the  towing  speed). 

ERROR  CHARACTERISTICS 

Individual  error  characteristics  suggest  how  they  can  be 
combined  into  an  overall  measurement  error  estimate.  How  often  an 
error  changes  in  value  is  one  such  characteristic.  Errors  are  clas¬ 
sified  as  "fast"  or  "slow  changing",  depending  on  how  frequently 
their  value  at  a  location  changes  significantly  during  a  typical 
measurement  time  (5  to  20  s) . 

The  contribution  of  "fast  changing"  errors  to  the  mean  measure¬ 
ment  error  gets  smaller  as  the  set  of  averaged  measurements  grows 
larger.  If  the  probability  density  function  of  the  "fast"  errors 
form  a  Guassian  distribution  about  a  mean  of  zero,  then  the  Central 
Limit  Theoreni  of  classical  statistics  applies.  The  standard  devia¬ 
tion  of  error  from  averaged  measurements  can  be  simply  related  to  the 
standard  deviation  of  individual  measurements  and  the  number  of  meas¬ 
urements  averaged,  as  follows: 


(4) 


^avg  ^ind  / 

where  N  ;  number  of  independent  measurements  averaged 

^avg*  Std.  Dev.  of  error  from  averaged  measurements 
^ind*  Std.  Dev.  of  individual  measurements 

"Slow  changing"  errors  have  one  value  for  all  individual  meas¬ 
urements  taken  at  a  given  location.  The  error  for  measurement  aver¬ 
ages  is  the  same  as  the  individual  measurement  errors  and  is 
independent  of  the  number  of  individual  measurements  averaged. 

Errors  are  also  classified  as  "position  dependent"  or  "position 
independent"  according  to  whether  their  magnitude  changes 
significantly  within  the  measured  flow  field.  Accuracy  effects  of 
position  dependent  errors  can  often  be  reduced  if  regions  of  large 
error  can  be  anticipated  or  detected.  LDV  operating  parameters  can 
be  changed  to  values  (not  applicable  to  the  flow  field  as  a  whole) 
which  will  enhance  accuracy  in  a  particularly  troublesome  region. 

Table  2  lists  possible  measurement  errors,  their  classification, 
and  the  inputs  that  determine  their  magnitude.  All  three  sets  of 
inputs  are  required  for  analysis  of  measurement  errors.  These  inputs 
are ; 

1.  Operating  Parameters 

2.  Flow  Field  Characteristics 

3.  LDV  Equipment  Characteristics 

INDIVIDUAL  PARTICLE  VELOCITY  ERRORS 

LDV  measurement  systems  really  measure  the  speed  of  small  parti¬ 
cles  traveling  with  a  fluid.  The  following  four  sources  of  error  can 
affect  the  accuracy  of  the  particle  velocity  component  measurements. 

Signal  Noise  Error 

Random  noise  in  the  scattered  light  signal  will  cause  velocity 
measurement  errors.  A  noise-free  signal  leaving  the  system 
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photomultiplier  has  nearly  a  sinusoidal  voltage  variation  (Figure  3). 
Real  signals  have  random  noise  (mean  value  of  zero)  added  from  vari¬ 
ous  light  and  electronic  sources.  The  LDV  system  signal  processors 
time  Mj  fringe  crossings  where  a  fringe  is  assumed  crossed  when  the 
voltage  signal  is  zero.  Errors  result  when  the  first  or  last  of  the 
timed  "zero-crossings"  is  advanced  or  retarded  by  the  signal  noise. 

Only  random  noise  signals  that  change  during  the  fringe  count 
will  cause  this  error.  So  this  error  must  be  considered  "fast  chang¬ 
ing".  Because  a  portion  of  the  random  noise  is  derived  from  noise  in 
the  received  scattered  light  this  error  source  is  "position 
dependent" .  Measurement  points  near  solid  surfaces,  for  example,  can 
be  expected  to  have  more  scattered  light  noise. 

Every  LDV  signal  burst  measured  has  its  own  SNR  or 
"signal- to-noise  ratio",  which  is  defined  here  to  be  the  ratio  of  the 
LDV  signal  amplitude  (peak)  to  the  standard  deviation  of  the  imposed 
random  noise.  Error  distributions  for  both  starting  and  for  stopping 
the  timing  of  fringe  crossings  are  equal  and  independent.  If  the 
LDV  signal  is  assumed  to  be  sinusoidal,  then  the  resulting  error  in 
particle  velocity  component  has  a  zero  mean  and  dimensionless 
standard  deviation  of: 

En  =  hf’sig*Sx/^o)  *  *  sin“^  d/SNR)  )  /  (2*tt*M^)|  (5) 

where  ^Sig*  measured  signal  frequency 

;  u,  V,  or  W  component  fringe  spacing 
Uq  :  towing  speed 

Note  that  increased  SNR  and  increased  fringe  counts  (M^:)  will 
reduce  this  error.  However,  simply  increasing  the  frequency  shift  to 
allow  counting  more  fringes  generally  does  not  reduce  the  error. 
There  is  a  canceling  effect  between  increased  values  possible  and 
the  increased  signal  frequency, 

Towing  tank  LDV  s igna 1- to-noi se  ratios  have  a  wide  distribution 
due  to  different  particle  sizes  and  trajectories.  They  also  differ 
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between  velocity  components  and  depend  on  counter  processor  adjust¬ 
ment  (Table  2).  Thus  this  formula  will  not  be  used  to  directly  esti¬ 
mate  signal  noise  error.  Rather,  experimental  particle  velocity 
component  data  (discussed  in  the  time  digitization  error  section) 
will  allow  a  calculation  of  typical  towing  tank  SNR  values  from 
measured  error  standard  deviations.  Table  3  lists  the  results. 

TABLE  3  -  DETERMINATION  OF  TYPICAL  COMPONENT  SNR  VALUES 


Variable 

Streamwise 

Component 

"U" 

Ve  1 1 i ca  1 
Component 

-•V  ” 

On-Ax  is 
Component 
"W" 

(E^2  ,  e^2)0.5 

0.0099 

0.0048 

0.0460 

Timing  Digitization 
Error:  E^ 

0.0002 

0.0002 

0.0011 

Signal  Noise 

Ector:  E^ 

0.0099 

0.0048 

0.0460 

^ig  *  ''  ^0 

1.000 

0.683 

4,78 

Typical  SNR  Values 
(eq.l  5) 

2.9 

4.0 

1.6 

Timing  Digitization  Error 

The  digital  data  transferred  from  the  LDV  signal  processors  to 
the  data  processing  computer  is  16-bits  long  with  a  12-bit  mantissa 
nci  a  -l-bit  exponent.  The  16-bit  word  represents  the  time  needed  for 
■  i  mrticle  to  pass  fringes  in  the  measurement  volume.  The  proces¬ 
sor  insures  (with  automatic  exponent  control)  that  the  12-bit 


mantissa  is  a  number  somewhere  between  2048  and  4096*  with  an  ac¬ 
curacy  of  +/-  1. 

The  digitization  error  changes  with  every  particle  measurement 
("fast  changing'*)  .  The  error  is  purely  random  and  independent  of 
fluid  mean  velocity  and  measurement  volume  location  ("position 
independent") . 

A  mantissa  error  that  is  equally  likely  to  be  +1,  0,  or  -1  is 
assumed.  The  mantissa  value  is  equally  likely  to  take  on  any  value 
between  2048  and  4096.  Then  the  digitization  error  in  particle 
velocity  components  has  a  mean  of  zero  and  dimensionless  standard 
deviation  of; 

4096 

=  (2/3)  *  (  Z  (l/I)  /  2048)  *  (Fg.  *S„/U_)  (6) 

1=2048  ° 

=  2.3*10"'^  *  (F_.  *s  /U  ) 

Sig  ^o' 

This  error  and  the  signal  noise  error  are  the  only  two  " fast 
chanq  inq"  errors  affecting  particle  velocity  component  measurements. 
Thus  in  a  steady  flow,  the  scatter  of  measured  particle  velocities 
about  the  true  mean  is  a  result  of  only  these  errors.  Typical  values 
for  the  standard  deviation  of  particle  velocity  components  appear  in 
Table  3.  These  were  taken  when  no  towed  model  was  present  to  disturb 
the  flow  or  introduce  noise  from  reflected  laser  beams.  The  towing 
speed  was  3  m/s  and  there  were  fixed  frequency  shifts  on  all 
components  (U:  0.0  MHz;  V:  0.5  MHz;  W;  0.5  MHz).  After  reducing 

the  error  due  to  digitization  effects  (insignificant),  typical 
component  SNR  can  be  calulated  using  Equation  (5).  These  also  appear 
in  Table  3.  Note  that  the  following  operating  parameter  choices  help 
increase  these  SNR  values  and  are  recommended  for  general  use; 

1.  Fringe  counts,  M^  >  4 

2.  Count  comparison  =  1% 


For  measured  frequencies,  Fg^g,  less  than  3.9  MHz 
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3. 

4. 


Signal  gain  holds  interburst  noise  to  within  +/-  25  mV 
Erroneous  data  eliminated  by  Chauvenet's  criterion 


Fringe  Spacing  Determination  Error 

Particle  speed  determination  require  an  accurate  value  for  the 
fringe  spacing.  For  a  properly  aligned  LDV  system*,  fringe  spacing 
is  uniform  and  depends  only  on  the  laser  light  wavelength  in  water 
and  the  laser  beam  crossing  angle.  The  laser  light  wavelengths  are 
known  to  within  +/-  .04%  (Appendix  B) .  Determination  of  the  beam 
crossing  angle  will  not  be  as  accurate. 

A  new  beam  crossing  angle  is  possible  only  when  the  LDV  optics 
are  realigned.  The  optics  should  not  be  realigned  during  the  course 
of  an  experiment.  Thus,  this  error  is  “slow  changing”  and  "position 
i ndependent . 

The  suggested  fringe  spacing  determination  technique  is  to  meas¬ 
ure  a  solid  surface  speed  (spinning  wheel)  underwater.  The  wheel  is 
spun  by  a  precisely  controlled  stepper  motor  (Figure  5a)  and  has  a 

rim  speed  of  0.2815  +/-  .0001  m/s.  This  error  estimate  results  from 

wheel  diameter  measurement  inaccuracies  (0.001  in.  or  0.05%  of  the 
diameter).  LDV  signal  noise  also  contributes  to  the  fringe  spacing 
error.  The  average  of  at  least  750  frequency  measurements  (with  no 
effective  frequency  shift)  are  recommended  for  the  vertical  and 
streamwise  spacing  determination. 

Sy  =  0.2815  (m/s)  /  Fy  (Hz)  (7a) 

=  0.2815  (m/s)  /  Fy  (Hz)  (7b) 

The  determination  of  on-axis  fringe  spacing  is  calculated  based  on 
the  streamwise  fringe  spacing  and  green  beam  wavelength  ( 
0.3859*10"^  meter  for  18°  C.;  sec  Appendix  B)  underwater. 


★ 

Fringe  spacing  can  change  across  the  measurement  volume 
if  beam  focal  points  and  crossing  points  do  not  coincide. 
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/  (2*sin(  sin  ^(^/2S^)/2  ))^  eq.#7c 

The  wheel-speed  error  and  exhibited  measured  frequency  repeat ibi 1 i ty 
(standard  deviation  of  0.06%  of  mean)  suggest  a  fringe-spacing- 
determination  error  with  mean  of  zero  and  standard  deviation  of: 

Streamwise  or  L 

Ef  =  7*10"^  *  U/Uq  (8a) 

Vertical  or  "V"  Component 

Ef  -  7*10"^  *  V/U^  (8b) 

On-A X  i_s  _qr  "  W ^qm p qn e n t : 

Ef  =  14*10"''  *  W/i'o  (8c) 

An  alternate  scheme  for  determining  fringe  spacings  is  to 
directly  measure  beam  crossing  angles  as  illustrated  in  Figure  5b. 
The  magnitude  of  error  in  the  angular  measurement  is  the  same  for 
green  and  blue  beam  pairs.  So,  fringe  spacing  determinations  for  the 
streamwise  direction  ( larger  X  value)  are  more  accurate  than  for  the 
blue  beams.  This  method  for  vertical  fringe  spacing  is  inadequate 
because  angular  det^' rm  i  na  t  i  on  errors  on  the  order  of  1%  are  possible. 
In  a  test,  the  streamwise  (or  green  beam)  fringe  spacing  determined 
in  this  way  w:.:-  only  0.2%  different  from  tne  spacing  found  from  the 
spinning  wheel.  In  the  same  test,  the  vertical  fringe  spacing 
determinations  liiffored  by  0.5%. 

Frequency  Sh.  i  f  t  Value  Error 

TSI,  Inc.  (the  LDV  system  niinu  facturer )  rates  their  effective 
frequency  sh!f*:s  in  be  .iceurate  to  within  r/-  0.01%.  The  frequency 
shifr  mus^  Oe  sub'i  acteo*  ‘  rom  the  measured  sugnal  frequency  before 
the  remainder  is  mu  I  t  iplied  by  t-n,.  fringe  spacing  to  obtain  the 
particle  speed. 

It  is  unc  M  t  lin  how  last  th'^  frequency  shift  error  changes.  It 
will  be  assumed  to  be  "slow  changing"  because  of  aging  electronic 
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components.  The  error  is  "position  independent”  because  it  is  purely 
a  result  of  electronic  component  operation. 

Assuming  a  uniform  error  distribution  between  +0.01%  and  -0.01%, 
the  frequency  shift  error  in  particle  velocity  components  has  a  mean 
of  zero  and  a  standard  deviation  of: 

E3  -  1.0001/1.73)  •  (Fsiq*S^/Uo)  (9) 

-  5.8*10-5  •  (Fsig.S,/U„) 


,-u* 


oJ 


r 


•  ^ 


Summary  of  Particle  Velocity  Errors 

The  four  particle  velocity  component  errors  of  this  section  are 
independent.  Three  of  the  errors  depend  similarly  on  signal 
frequency  measured  and  the  component  fringe  spacing  term 
(E’g j^g  +  Sjj/UQ)  .  This  is  because  these  errors  directly  affect  the 
measured  signal  frequency,  and  not  the  frequency  directly 
proportional  to  particle  velocity,  (F„.  -ps) . 

w  X  ^  ' 

Two  of  the  errors  are  "fast  changing"  and  will  be  reduced  by 
taking  the  mean  of  multiple  samples  of  the  particle  velocity.  The 
two  "slow  changing"  errors  will  not  be  reduced  in  multiple  sample  me¬ 
ans.  In  most  situations  and  E^  will  far  outweigh  the  other  two 
errors.  Table  4  conveniently  displays  this  information. 

MEAN  FLUID  VELOCITY  ERRORS 

The  towing  tank  LDV  system  was  intended  to  obtain  the  3 
component,  mean  water  velocity  in  the  flow  field  around  ship  models. 
This  section  describes  additional  errors  introduced  when  turning  the 
particle  velocity  measurements  of  the  previous  section  into  mean  wa¬ 
ter  velocity  estimates. 

Particle  "Lag"  Error 

This  error  represents  the  difference  between  particle  velocity 
and  the  local  water  velocity.  The  ability  of  a  particle  to  "follow" 
precisely  the  local  flow  field  is  a  function  of  the  particle's  size 
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TABLE  4 


PARTICLE  VELOCITY  COMPONENT  ERROR  FORMULAS 


Character 

Fast  Slow 

Streamwi se 

Component  Velocity  Errors 

Ve  r  1 1 ca  1 

On- Ax  i  s 

X 

0.  132*F/(Ug»Mf) 

0.230*F/(Uo*Mf) 

4.42*F/(Uo 

X 

.0004*F/ (Ug*Mf ) 

.0010*F/ (Uo*Mf > 

.0066*F/  (U, 

X 

.0007*1U/U,j| 

.0007* |V/Uo 1 

.0014* |W/U, 

X 

.0001»F/Uo 

.  000  3  *F/Uo 

.0017*F/Ug 

KEY  : 

F  :  coTiponent  signal  frequency  in  MHz 
U^:  towinq  speed  in  m/sec 

signal  processor  fringe  count 


and  v/eight  and  of  water  accelerations  near  the  measurement  point. 
Figure  6  shows  the  density  and  size  distribution  of  the  seed  parti¬ 
cles  initially  chosen  for  towing  tank  LDV  measurements. 

The  amount  of  particle  lag,  at  a  given  measurement  location,  can 
change  as  fast  as  the  velocity  field  ("fast  changing").  However,  in 
the  presence  of  some  persistant  fluid  flow  structures  (e.g.,  a  trail¬ 
ing  vortex  from  a  lifting  surface,  the  expected  value  of  the  error  is 
not  zero  ("slow  changing") .  Particle  lag  errors  are  " pos i t i on 
dependent"  because  they  result  from  the  local  fluctuating  and/or  mean 
water  velocity. 

Several  analyses  concerning  particle  lag  are  done  in  Appendix  C. 
rh<_'  first  analysis  considers  the  particle's  local  velocity  to  have  a 
steidy  mean  and  a  sinusoidal  variation  with  time. 


Particles  tend  to 


silver  Coated  Class  Spheres  (5  ■  2.8  gm/cc) 


Figure  6  Diameter  Distribution  of  Seed  Particles 


follow  this  fluctuation  due  to  viscous  forces*.  Stoke' s  flow  about  a 
spherical  particle  is  assumed.  The  particle  follows  the  flow  with  a 
sinusoidally  varying  velocity  that  differs  in  amplitude  and  phase 
from  the  water  velocity.  The  difference  is  a  function  of  the  parti¬ 


cle  size  and  density.  Figure  7  shows 
an  assumed  fluctuation  magnitude  and  a 


the  amplitude  error  for 

fluctuation  frequency  of  1 


KHz.  It  is  based  on  the  following  equation  derived  in  Appendix  C. 


3.88*10 


■2  *p  2  * 

Vel 


(Uo  *y  ^ 


where  ^Vel*  frequency  of  velocity  fluctuations 
SGp  :  particle  specific  gravity 
Dp  :  particle  diameter 

(T^  :  std.  dev.  of  velocity  fluctuation 

V  ;  kinematic  viscosity  of  water 


Low  Reynold's  numbers  result  from  the  particle  size. 
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Assumptions : 

SG-  •  2.8 


1000  Hi 
3.0  m/s 


TOWINO  TA,NK  SEED 
Si'zE  DISTRIBUTION 


Particle  Diameter  (i/m) 


Figure  7;  Estimates  of  Particle  Lag  Error 


A  second  analysis  (Appendix  C)  examines  particle  lag  in  vortical 
flows  where  the  difference  in  fluid  and  particle  velocities  is  found 
as  a  function  of  expected  particle  seed  accelerations.  The  accelera¬ 
tion  due  to  large  scale  eddies  (length  scale  of  wake  boundary  layer 
thickness  and  velocity  scale  of  the  towing  speed)  and  the  smallest 
detectable  eddies  (length  scale  of  LDV  measurement  volume  are  both 


cons ide  red 


It  is  shown  that  the  smaller  scale  eddies  in  this 


analysis  cause  the  most  particle  lag.  Figure  7  contains  curves  show¬ 
ing  the  results  of  the  lag  error  calculation  for  certain  assumptions, 


It  IS  based  on; 


\2  -  («o  *  SGp*Dp2)  /  (X8*V.  (t2*2*D^„)»-“) 


where  S  :  boundary  layer  thickness 

No  estimate  of  error  is  made  for  persistant,  large  vortices. 
These  can  be  evaluated  only  after  some  velocity  measurements  of  the 
particular  flow  structure  have  been  taken.  Also,  note  that  because 
of  this  centrifugal  effect,  the  inner  radii  of  vortices  will  have 
fewer  and  smaller  particles.  This  could  be  a  source  of  "bias"  in  ad¬ 
dition  to  those  listed  in  Table  2. 

A  generally  applicable  estimate  of  particle  lag  error  magnitude 
can  be  made  by  arbitrarily  combining  the  two  separate  lag  error 
estimates : 


=  (  (3.9*10“'^*  ^T/Uq)  2  +  (3.4*10‘'^*Uq)  2  )0-5 


where  U^,  (T^ :  velocities  in  m/s 
Eyei  =  1000  Hz 
Dp  =  6 . 0  urn 

SGp  .  2.8 

i  =  0.15  m 


Velocity  Fluctuation  Error 

This  error  occurs  because  the  LDV  averaging  time,  T  „  is 

^  ^  avg 

shorter  than  the  time  needed  for  significant  velocity  fluctuations. 
Velocity  measurements  differ  from  the  desired  mean  value  because  of 
turbulence  and  blade  rate  fluctuations  near  ship  propellers.  This 
error  is  "fast  changing"  because  velocity  fluctuations  are  expected 
to  occur  much  faster  than  the  measurement  period  (2  -  10  seconds). 


Both  the  frequency  and  magnitude  of  velocity  fluctuations  will  change 
with  measurement  position  ("position  dependent”) . 

An  estimate  of  the  error  standard  deviation  must  come  from 
knowledge  of  the  local  velocity  fluctuations. 

^  =  <^x/Uo  (13) 

where  (T^ :  standard  deviation  of  fluctuating  velocity  component 


Bias  Errors 

The  towing  tank  LDV,  as  a  measurement  technique,  samples  the  wa¬ 
ter  velocity  only  when  a  particle  passes  through  the  measurement 
volume  and  crosses  fringes.  This  complicates  the  determination  of 
mean  velocities  because  some  velocity  magnitudes  and  directions  will 
tend  to  offer  up  measureable  particles  at  a  greater  rate  than  others. 
The  result  is  that  in  a  fluctuating  velocity  field,  a  simple  average 
of  particle  measurements  is  biased  from  the  true  mean  velocity  toward 
these  more  "favorable"  velocity  magnitudes  and  directions. 

Because  bias  errors  are  errors  in  mean  velocity  estimates,  they 
are  "slow  changing".  Their  dependence  on  the  magnitude  of  velocity 
fluctuations  and  the  distribution  of  particles  makes  them  "pos i t i on 
dependent" . 

Much  bias  is  eliminated  if  the  particles  are  distributed 
uniformly  in  the  measured  fluid.  The  degree  to  which  particle 
uniformity  is  achieved  in  the  towing  basin  under  standard  operating 
procedures  (i.e.,  seeding  only  on  return  passes  and  15  minute  towing 
cycles)  is  not  quantitatively  known.  Seed  is  sprayed  uniformly  over 
a  0.5  m.  wide  strip  centered  in  the  range  of  measurement  locations 
and  reaching  the  length  of  the  basin.  Seed  is  sprayed  nearly  every 
return  pass,  so  that  seed  falling  from  the  measurement  region  is 
gradually  replaced  and  not  introduced  in  large  concentrations. 
Obst' rvat  1  ons  of  seed  across  the  entire  basin  width  have  been  made 
afrer  only  4  seeding  cycles.  It  is  believed  that  these  procedures 
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and  observations  ensure  uniform  seed  distribution  vertically  and 
streamwise.  There  probably  is  a  measureable  seed  concentration  from 
side  to  side,  but  not  a  significant  one  over  the  central  meter  wide 
volume  of  water  that  can  conceivably  pass  through  the  LDV  measurement 
volume.  Aquisition  data  rates  (away  from  solid  surfaces)  have  not 
shown  any  noticeable  dependence  on  measurement  location. 

Even  if  particle  distribution  biases  are  insignificant,  poten¬ 
tially  important  biases  remain.  Dimotakis^^  derives  an  expression 
for  the  mean  rate  that  particles  moving  with  a  given  velocity  will 
produce  an  LDV  measurement.  Appendix  D  extends  Dimotakis'  work  to 
the  3  component,  frequency  shifted  LDV  towing  tank  system.  That  mean 
measurement  rate,  calculated  by  knowing  instantaneously  all  three 
components  of  velocity,  can  be  used  to  weight  individual  particle 
velocity  measurements.  The  weighting  allows  an  unbiased  estimate  of 
the  mean  fluid  velocity  (if  all  occurring  velocities  are  measureable 
to  some  extent) . 

The  calculated  mean  measurement  rate  of  Reference  11  is 

conveniently  split  into  the  product  of  two  functions  of  the  particle 
velocity.  The  first  function  G^(u)  is  the  mean  number  of  particles 
passing  through  the  measurement  volume  per  unit  time.  The  second 
function,  Gj(U),  is  the  fraction  of  the  arriving  particles  that 
produce  a  measureable  signal  (i.e.,  those  that  cross  a  minimum  of 
fringes) .  Velocity  measurements  normalized  by  are  free  of  what  is 
usually  called  "velocity  bias".  Velocity  measurements  normalized  by 
the  second  function,  are  free  of  what  has  been  called  "fringe 

bias"  . 

(  U.  /(Gy  (U  )  *  Gf  (U  )  )  ) 

U  .  -  (14)  I 

(l/(Gy  (U  )  *  Gf  (U  )  )  )  j 

I 

I 

I 

A  third  type  of  bias  results  if  frequency  filters  (on  the  i 

counter  signal  processors)  attenuate  some  valid  LDV  scattered  light 
signals  more  than  others.  The  counter  processor  will  record 
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particular  set  of  mean  velocities  and  a  range  of  velocity  fluctuation 
magnitudes.  The  velocity  fluctuation  for  each  component  is  assumed 
to  be  a  sinusoidal  variation  on  top  of  the  mean  velocities  listed  in 
the  figure.  Random  phase  means  that  the  phase  difference  between  the 
velocity  fluctuations  of  the  streamwise  component  ("U")  and  the 
transverse  components  ("V"  and  "W”)  was  varied  from  0°  to  360°.  The 
average  error  was  graphed.  When  there  was  no  phase  difference 
between  components,  the  errors  for  all  components  are  close  to  the 
random  phase  error  for  the  streamwise  component  ("U"). 

Velocity  bias  can  be  eliminated  by  weighted  data  averages  as 
shown  in  Appendix  D  and  Reference  11.  The  weighting  function,  is 
the  average  particle  arrival  rate  for  the  measured  velocity. 
However,  to  maximi2e  the  data  acquisition  rate,  the  towing  tank  LDV 
system  is  not  operated  to  acquire  all  three  velocity  components 
simultaneously.  Thus  it  is  not  possible  to  calculate  which 
requires  the  value  of  all  components  at  the  time  of  measurement. 

However,  proper  selection  of  counter  processor  parameters  allows 
another  estimator  of  average  particle  arrival  rate.  The  rate  at 
which  particles,  with  a  particular  velocity,  move  through  the  meas¬ 
urement  volume  is  equal  to: 

Particle  Rate  =  n  *  Vol  /  Tg(u)  (15) 

where  n  :  particles  per  unit  volume 

Vol  :  volume  of  measurement  region 

Tg(u)  :  average  residence  time  of  particles  with  velocity  U 
in  the  measurement  volume 

oince  the  particle  density  and  measurement  region  volume  are 
constants,  velocity  bias  can  also  be  eliminated  by  using  a  weighting 
function  of  "l/^g(n)".  The  towing  tank  LDV  processors  can  record  the 
particle  residence  time  Tg^  for  each  measurement  if  operated  in 
TBC  or  TBM  modes.  TSI  Inc.  suggests  its  use  as  a  weighting  function 
in  the  elimination  of  velocity  bias.^® 
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similarly  for  "V"  and  "W” 

The  processors  detect  the  particle  residence  time  by  counting 
fringe  crossings.  Random  noise  spikes  can  easily  cause  the  processor 
to  prematurely  think  that  a  particle  has  exited  the  measurement 
volume.  Residence  time  measurements  or  total  fringe  counts  will  show 
more  low  fringe  counts  if  this  effect  is  important.  Figure  9  shows 
the  distribution  of  measured  total  fringe  counts  (streamwise 
component)  and  the  distribution  expected  from  gaussian  laser  beams, 
the  seed  particle  distribution  (Figure  6  and  scattered  light  inten¬ 
sity  proportional  to  particle  diameter  squared)  and  measurement 
volume  geometry.  The  large  number  of  extra  low  fringe  counts  is 
obvious.  It  appears  that  as  many  as  half  of  the  total  fringe  counts 
are  mistakenly  low.  The  velocity  bias  correction  essentially  has  no 
effect  on  these  measurements,  so  a  reasonable  assumption  (for  Figure 
9  data)  is  that  as  much  as  half  of  the  velocity  bias  is  not  removed 
by  the  correction. 

It  is  shown  in  Appendix  E  that  Tg^(U)  has  a  standard  deviation 
about  Tg(U)  of  as  much  as  0.35T^g.  Appendix  E  also  shows  that  this 
results  in  a  standard  deviation  of  weighted  velocity  measurements 
about  U  that  is  as  much  as  6%  larger  than  the  true  ("^3)  weighted 
velocity.  So,  velocity  bias  correction  by  the  Tg^  weighting  method 
increases  velocity  fluctuation  standard  deviations  by  a  factor  of 
1.06  and  eliminates  one  half  (arbitrarily  chosen  based  on  Figure  9) 
of  the  errors  graphed  in  Figure  8. 

=  1.06  *  /  []o  ill) 

=  1.06  * 
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ge  Crossings  in  Measurement  Volume  (Streamwise  Component) 


I 

Fringe  Bias  Estimation  In  a  given  measurement  situation,  fringe  bias 
occurs  because  certain  velocities  are  more  likely  to  result  in  valid 
I  processor  measurements  than  other  velocities,  A  simple  mean  of  valid 

'  measurements  is  biased  toward  the  more  "measureable"  velocities.  The 

measureabi 1 i ty  of  any  velocity  is  reflected  by  the  ratio 
;  where  is  the  number  of  fringe  crossings  for  a  particle  with 

j  velocity  and  passing  through  the  measurement  volume  center. 

'  Apor ox imate ly : 


(U) 

(V) 
Nf  (W) 


D^V*(FS  +  u/s^) 
Dj^^MFS  +  V/Sy) 
D^ivMFS  +  W/S^) 


/  (U^  +  v2)0*  5 

/  (U^  +  ^  ^ 

/  (U^  +  )  0  •  5 


1 

whe  re 


D 

mv 

FS 
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diameter  of  component's  measurement  volume 
effective  frequency  shift  of  component 
fringe  spacing  of  subscripted  component 


(18a) 

(18b) 

(18c) 


Appendix  D  formulas  for  fringe  bias  require  measurement  volume 
dimensions.  These  can  be  arbitrarily  set  by  calculated  points 

(Figure  4),  However,  when  correctly  set,  the  boundaries  are  points 
wiiere  the  particle  scattered  light  just  produces  a  signal  above 
processor  thresholds.  The  boundaries  are  therefore  really  functions 
of  particle  size  and  processor  signal  gain  (as  controlled  by  signal 
noise  amplitude). 

An  experimental  determination  of  the  streamwise  measurement 
V  )  I umc  diameter  was  made  by  recording  the  total  number  of  detected 
fringe  crossings  for  an  essentially  streamwise  towing  tank  velocity 
prob('  strut  was  towed  without  a  model).  Figure  9  displays  a 
hi.s*ogriT'  of  the  results.  Signals  with  less  than  8  fringe  crossings 
r'  nri*-  recorded.  Also,  signals  with  noise  that  interrupted  the  to- 
il  fringi-  cnun*-  appear  as  low  fringe  counts.  Most  importantly 
'nowo-t-r  ,  an  upper  limit  to  the  number  of  fringe  crossings  is 
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apparent . 

Boundaries  of  D^_2  would  result  in  a  maximum  of  64  fringe  cross¬ 
ings  =  0.11  mm;  Sy  =  1.72  um)  .  The  figure  shows  as  many  as  113 
fringe  crossings.  However,  one  might  choose  71  fringe  crossings  as 
the  maximum  for  an  average  sized  particle  (90%  of  the  data  fall  at  or 
below  71  fringe  crossings). 

The  experimental  data  indicates  a  "detectable"  measurement 
volume  11%  larger  than  that  calculated  based  of  2*  Though  the 
experimental  data  is  only  for  the  streamwise  component,  it  is  assumed 
that  the  measurement  volume  dimensions  and  Lrnv^ 

components  are  1.11  times  greater  than  the  calculated  2  dimensions 
of  f ig ure  4  . 

Figure  10  displays  calculated  fringe  bias  errors  as  a  function 
of  velocity  fluctuation  magnitude.  Only  the  fringe  bias  errors  of 
tile  vertical  and  on-axis  components  are  large  enough  to  appear  on  the 
graph.  For  the  same  frequency  shift,  the  vertical  error  is  much  more 
important  than  the  on-axis  component  error.  This  is  generally  true 
for  the  velocities  and  frequency  shifts  found  on  towing  tank  LDV 
experiments.  The  velocity  means  and  fluctuation  characteristics  are 
the  same  as  in  figure  8  for  velocity  bias.  For  fringe  bias,  the 
choice  of  frequency  shift  (affecting  N^:)  and  Mj-  values  are  important. 
This  is  because  fringe  bias  errors  vary  directly  with  the  mean  and 
the  standard  deviation  of  the  ratio  that  occur  while  measuring 
a  component  mean. 

No  direct  formula  can  be  given  for  unco  r  rt'^c  ted  fringe  bias  er¬ 
rors.  However,  with  some  assumptions  about  the  nature  of  the 
r't'locity  fluctuations,  t'no  calculations  finure  10  can  be 
r>-’pea  t.-.'d  (Appendix  D).  In  this  way  an  error  ('stimatr'  can  be  made. 

^1  tcM.'  Bias  Estimation  It  is  very  possible  tti  ,  i  lie  desire  to  reduce 
si'jnal  noise  and  increase  data  rates  in  th'-  lowing  "ink  LDV  system 
will  rule  against  setting  a  large,  saf(’  (with  r-'spect  to  filter  bias) 
ringe  of  unfiltered  frequencies  on  the  ciounti'i  processor.  Rather, 
the  low  pass  and  high  pass  filters  will  be  e'C  e,  -Lose  together  as 


still  allows  measurement  of  all  velocities  that  occur.  Since  high 
frequency  noise  is  the  most  disruptive  of  counter  processor  opera¬ 
tions,  the  low  pass  filter  is  particularly  important.  Appendix  F 
estimates  the  amount  of  bias  developed  as  the  signal  frequency 
approaches  the  low  pass  filter  frequency.  Counter  processor  specifi¬ 
cations  rate  these  filters  with  a  minimum  of  30  dB  per  octave  rol¬ 
loff. 

Figure  11  represents  the  results  of  Appendix  F  under  the  follow¬ 
ing  assummptions : 

1.  low  pass  filter  rolloff  of  30  dB  per  octave  starting  at 

the  filter  frequency,  F„.,^ 

^  ^  '  Filter* 


s 


les 


a  mean 


Sig 


and  standard  devia¬ 


tion 


Sig 


3.  the  signal  frequencies  have  a  normal  distribution. 

4.  frequency  measureabi 1 i ty  is  directly  proportional  to  the 
precent  of  the  signal  voltage  passed  by  the  low  pass  filter. 

Note  that  the  error  given  in  the  figure  is  a  percentage  of  F  . 

O  1  ^ 

To  get  the  filter  bias  velocity  error  as  a  percent  of  the  towing 
speed,  this  number  must  be  multiplied  by: 


K  =  F^.  *  e 

Sig 


/  u, 


(19) 


Summary  of  Mean  Fluid  Velocity  Errors 

All  the  errors  in  this  section  are  assumed  to  be  independent  of 
each  other  and  all  previously  defined  errors.  Formulas  are  developed 
for  all  the  errors  except  for  the  three  bias  errors.  These  three  can 
also  be  estimated  from  operating  parameters  and  flow  field 
characteristics.  However,  the  nature  of  these  calculations  requires 
a  computer  program  to  do  the  long  summations  and  to  tailor  the  many 
flow  field  assumptions  to  the  user's  requirements.  With  the  avail¬ 
able  velocity  bias  correction,  the  bias  is  "half"  removed  but  a 
modified  velocity  fluctuation  error,  is  required. 

Two  of  the  errors  are  "fast  changing"  and  will  be  reduced  by 
taking  the  mean  of  multiple  samples  of  the  particle  velocity.  This 
is  d  good  assumption  for  the  type  of  particle  lag  error  formulated. 
However,  particle  lag  error  due  to  some  persistant  fluid  flow 
structures  is  "slow  changing"  and  will  not  be  accounted  for.  The 
user  must  be  on  the  look  out  for  these  situations  and  estimate  his 
own  particle  accelerations  and  lag  errors  (using  Appendix  C  formu- 
l  a  s  )  . 

The  uncorrected  bias  errors  are  "slow  changing"  errors.  They 
will  not  be  reduced  in  multiple  sample  means. 

Tn  most  situations  (if  frequ(-'ncy  shifts  are  properly  chosen  and 


velocity  bias  correction  is  done) ,  will  be  the  only  important  er¬ 
ror  of  this  group.  Table  5  conveniently  displays  this  information. 


TABLE  5  -  MEAN  WATER  VELOCITY  ERROR  FORMULAS 


Character 


Component  Velocity  Errors 


Fast  Slow 


Streamwise 


Vertical 


On-Ax  IS 


3.4*10'  ‘’'((‘^x/Oo*^  ^ 


/  Uo 


"“v  /  Oc 


/  “o 


"■VB- 
Uncor . 


E'  =  1.06*E„  Ei  =  1.06*E„  Ei  =  1.06*E^ 
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:  towing  speed 

'^x  '  *^J  '  ^V' standard  deviation  of  velocity  fluctuation 

C.P.  :  computer  program  for  estimation 

‘’"sig  :  standard  deviation  of  signal  frequency 


FINALIZED  WATER  VELOCITY  DATA  ERRORS 


Mean  water  velocity  components  already  analyzed  are  not  the  end 
result  desired  from  typical  wake  surveys.  Ideally  the  user  seeks,  at 
known  positions,  non-dimensionalized,  velocity  components  that  are 
orthogonally  aligned  to  the  ship  model  and  undisturbed  by  the  LDV 
strut.  Velocities  are  usually  non-dimensionalized  by  an 
independently  measured  free  stream  velocity  magnitude  (towing  speed). 
Component  directions  are  not  exactly  orthogonal  due  to  optical  align¬ 
ment  limitations.  Velocity  disturbances  created  by  the  LDV  strut  and 
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positioning  inaccuracies  cause  measured  velocities  to  be  different 
from  those  sought.  These  four  error  sources  and  available  correction 
techniques  are  treated  in  this  section. 

Carriage  Speed  Errors 

This  error  is  the  result  of  towing  carriage  speed  variation.  It 
is  not  an  error  in  fluid  velocity  measurement  Rather,  the  error  is 
introduced  when  the  measured  data  is  non-d imensional i zed . 

Many  carriage  passes  will  be  required  to  take  an  entire  wake 
survey.  In  order  to  use  data  taken  on  different  carriage  passes  or 
even  at  different  times  on  the  same  carriage  pass,  the  measurements 
must  all  be  non-d imens i ona 1 i zed  by  the  carriage  speed  at  the  time  of 
measurement . 

The  LDV  measurement  procedure  records  an  average  carriage  speed 
for  each  pass.  It  uses  that  speed  to  non-dimensional ize  every  meas¬ 
urement  made  during  the  pass.  Thus,  this  error  is  "pos i t i on 
i ndependent”  and  "slow  changing”.  Speed  variations  during  the  same 
pass  can,  with  care,  be  held  to  within  +/-  .1%  of  the  mean  speed. 
This  is  better  than  the  carriage  speed  variation  between  different 
passes  that  can  be  +/-  .3%  or  more  when  equal  speed  passes  are 
sought.  Assuming  a  normal  distribution  of  towing  carriage  speed 
(during  a  pass)  with  a  mean  of  and  standard  deviation  of 
5*10  carriage  speed  error  has  a  mean  of  zero  and  a  standard 
deviation  of: 


5*10' 


Carriage  speed  measurement  errors  and  flow  field  changes  (from  varia¬ 
tion  of  model  ship  Froude  number)  have  been  neglected  in  this  error 
.?  s  t  i  m  a  t  i  o  n  . 

Flow  Disturbance  Errors 

The  configuration  of  the  towing  tank  LDV  system  means  that  the 
partial  Ly  submerged  probe  strut  (figure  2)  will  always  cause  some 


disturbance  of  the  flow  field.  The  strut's  location  and  basic  shape 
are  mostly  set  by  the  optical  requirements  of  the  laser  beams  and 

scattered  light  that  pass  through  it.  In  order  to  provide  for 
rigidity  and  laser  light  passage,  the  underwater  portion  of  the  strut 
is  basically  a  hollow  aluminum  column;  0.19  m.  wide  by  0.095 
m.  thick  by  as  much  as  0.90  m.  tall.  The  lowest  portion  (0.18  m.)  of 
the  strut  must  accommodate  the  laser  beam  focusing  lens  and  is  0.122 
m.  thick. 

Plastic  fairings  are  attached  to  the  forward  and  rear  faces  of 
the  strut  to: 

1.  reduce  separation  for  flow  angles  of  -5°  to  5°; 

2.  reduce  strut  drag; 

3.  reduce  flow  disturbance  at  the  LDV  measurement  volume; 

The  faired  shape  appears  in  figure  12.  The  asymmetrical  shape  of  the 
box  fairing  causes  less  flow  disturbance  on  the  flat  side  (side  of 
LDV  measurement  volume)  than  the  more  curved  side.  However,  it  also 
produces  a  sideways  force  or  lift  that  may  cause  strut  deflection  and 
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Figure  12:  LDV  Strut  Dimonsions 


stability  problams  above  the  maximum  checkout  speed  of  4.0  m/s. 

Several  experimental  tests  were  run  to  quantify  the  flow 
disturbance  of  the  strut.  Basically  they  each  consisted  of  putting 
another  velocity  sensing  instrument  (pitot  static  or  5-hole  pitot 
tube)  at  the  LDV  measurement  volume  location.  On  a  single  carriage 
pass,  the  pitot  tube  measurement  would  be  recorded,  the  LDV  strut 
would  be  traversed  out  of  the  water  (or  0.6  m  to  the  side),  and  a 
second  pitot  tube  measurement  would  be  recorded.  The  difference 
between  the  first  and  second  pitot  tube  measurements  is  reported  here 
as  the  LDV  probe  strut  flow  disturbance. 

The  first  series  of  measurements  was  done  in  the  free  stream 
with  no  ship  model.  Potential  flow  calculations  indicated  that  the 
streamwise  velocity  component  would  be  disturbed  the  most.  A  pitot 
static  tube  that  sensed  only  that  velocity  component  was  used.  The 
results  appear  in  figure  13  for  a  range  of  strut  depths.  is  the 
pitot  tube  measurement  of  the  disturbed  flow  with  the  LDV  strut  in 
place.  All  disturbance  values  are  less  than  0.4%  of  the  tow  speed. 
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velocity  for  both  the  disturbed  and  undisturbed  cases.  Considering 
the  data  scatter  (especially  near  the  angular  location  of  180°) ,  the 
undisturbed  and  disturbed  data  sets  are  essentially  the  same  cr  at 
worst  different  by  about  .01*0^.  The  very  wide  model  (an  auxiliary 
oiler)  and  the  strong  (U  <  0.2*0^),  unstable  wake  field  at  the  180° 
location  were  purposely  chosen  as  a  severe  test  of  the  potential  for 
strut  flow  disturbance.  Most  wake  survey  positions  will  have  flow 
disturbances  closer  in  magnitude  to  those  found  in  the  free  stream 
tests  (  figure  13 ) . 

Typical  error  magnitudes  are  illustrated  in  figures  13  -  15.  It 
is  suggested  that  the  strut  disturbance  correction  (discussed  in  the 
next  section)  always  be  used.  This  correction,  based  on  free  stream 
velocity  measurements,  eliminates  free  stream  flow  disturbances  for 
all  velocity  components.  For  error  analysis,  the  flow  disturbance 
error  will  then  be  assumed  to  be  zero.  In  reality  though,  the  LDV 
user  must  still  be  wary  of  situations  like  measuring  "unstable", 
"severe  wake"  flows  that  are  we^l  under  "wide  models". 


E.  =  2*10 


(uncorrected) 

(corrected) 


Velocity  Component  Directional  Errors 

Despite  careful  alignment  procedures,  the  three  measured 
velocity  components  do  not  come  out  perfectly  perpendicular  and 
aligned  to  the  ship  model.  Errors  of  as  much  as  .5°  are  possible  for 
the  vertical  and  streamwise  components.  The  on-axis  component  can  be 
as  much  as  1.5°  misaligned.  These  are  "slow  changing"  and  "pos i t ion 
i ndependent"  errors.  Special  caution  should  be  observed  with  the 
on-axis  component  because  even  the  slightest  change  in  optical  align¬ 
ment  can  have  a  large  effect  on  this  error.  The  correction  scheme 
described  in  this  section  should  be  done  after  any  optical  change  to 
the  LDV  system.  As  a  further  precaution  it  should  be  done  every  week 
of  experimentation,  even  if  there  is  no  known  optical  change. 


ir 


As  long  as  the  measured  velocity  components  are  not  coplanar, 
the  desired  orthogonal  velocity  components  can  be  calculated  if  the 

’""easured  directions  are  known.  Alignment  procedures  make  all  direc¬ 
tional  errors  zero  to  the  best  of  their  limited  precision.  Further 
(but  not  complete)  correction  is  possible  from  data  taken  on  calibra¬ 
tion  runs  with  a  constant  carriage  speed  and  no  model. 

For  a  perfectly  aligned  and  otherwise  accurate  LDV  system,  the 
measured  velocity  components  on  calibration  runs  should  be; 


U  =  *  cos  (0) 

V  =  Uq  *  sin  (0) 

W  =  0.0 

where  0  =  vertical  traverse  angle 


(23b) 

(23c) 


The  difference  of  measurements  from  these  values  is  the  result  of  a 
number  of  errors: 

1.  transverse  components  not  perpendicular  to  streamwise 
direction  (E^) 

2.  carriage  speed  determination  (E^) 

3.  streamwise  fringe  spacing  determination  (E^) 

4.  strut  flow  interference  (Ej) 

5.  frequency  shift  -  signal  mixing  accuracy  (Eg) 

Of  the  five  errors  listed,  only  E^  should  not  be  proportional  to 
the  towing  speed  (should  change  randomly  with  channel  frequency 
shift)  .  There  is  no  way  and  little  need  to  separate  the  proportional 
errors.  If  calibration  runs  are  done  for  a  range  of  tow  speeds 
(expected  U  velocity  component  range),  corrections  proportional  to 
tow  speed  and  constant  corrections  can  be  separated  resulting  in  six 
correction  factors: 


U  *  cos(0)  = 


Au*U,  +  Bu 


Uq  *  sin(0)  =  +  Ay*U  + 


(24b) 


■/  f  ^ 


(24c) 


•  V-i-' 


Wj  +  A^J*U  + 


where 


Ur,Vj.,Wi.: 
®U » By , : 


raw,  uncorrected  velocity  components 
proportional  correction  factors 


constant  correction  factors 


In  general  the  B  correction  factors  should  be  small  enough  to 
neglect.  Doing  only  one  tow  speed  and  assuming  B^  =  0  is  just  about 
as  good  as  calculating  the  more  involved  two  parameter  correction. 
If  parameters  are  derived  from  many  points  on  several  runs,  the  ef¬ 


fect  of  tow  speed  measurement  errors  (mean  value  of  =  O)  will  be 


reduced  through  averaging.  If  parameters  are  derived  from  points  at 
various  depths  (in  the  range  used  during  the  experiment)  then  an 
average  flow  disturbance  will  be  corrected.  The  experimental  data 
correction  formulae  then  are: 


+  Bu 


(25a) 


V  =  Vr  +  Av*U  +  B^ 


(25b) 


W  =  +  Am*U  +  Bi 


(25c) 


This  correction  does  nothing  about  streamwise  component  direc¬ 
tion  or  vertical  component  and  on-axis  component  directions  within 
the  transverse  plane.  Sizable  alignment  errors  of  this  type  are 
still  possible  (0.5°  for  vertical  or  streamwise  and  1.5°  for 
on-axis).  The  estimated  directional  alignment  errors  that  follow  are 
based  on  the  alignments  left  uncorrected  and  on  the  repeatability  of 
the  calculated  correction  factors.  Also,  note  that  full  alignment 
correction,  means  that  previous  error  standard  deviation  magnitudes 
should  be  assumed  zero  (E^,  -  streamwise  only,  E^ ,  and  Eg). 
Uncorrected  and  corrected  alignment  directional  errors  should  have  a 
mean  of  zero  (over  many  optical  system  setups)  and  a  standard  devia¬ 


tion  of : 


•  •  •  I 
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No  Correction  Factors  Used 


£3(0)  =  ((.005*V)^ 

+ 

(.005*W)2  )0*5  / 

(26) 

EgCV)  =  ((.005*U)^ 

+ 

(.005*W)2  )0-5  / 

E3(W)  =  ((.013*U)^ 

+ 

(.013*V)2  /  U^ 

Correction  Factors  Used 

Ea(U)  =  ((.005*V)^ 

+ 

(.005*W)2  +  {.001*Uq) 2)0.5  / 

(27) 

EgiV)  =  ({.005*W)^ 

+ 

(.001*Uq)2)0.5  / 

E3(W)  =  ((.013*V)2 

+ 

(.002*U^)2)0.5  / 

Traverse  Positioning  Errors 

In  a  non-uniform  velocity  field,  positioning  errors  result  in 
measured  velocities  that  are  different  from  the  velocity  at  the  in¬ 
tended  position.  The  velocity  error  depends  on  both  the  magnitude  of 
the  positioning  error  and  the  magnitude  of  the  local  velocity 
gradient.  Five  sources  of  positioning  error  are  discussed  below. 

Minor _ Positioning  Errors.  The  SONY  Magnescale  encoders  incorporated 

into  the  LDV  traverse  system^^,  have  an  accuracy  of  +/-  .003  mm. 
Positioning  is  automatically  checked  by  these  encoders  before  and 
after  each  carriage  pass. 

During  the  carriage  pass  the  traverse  system  stepping  motors  are 
relied  on  to  make  the  correct  position  changes  according  to  a 
calibrated  ratio  of  motor  steps  to  position  change.  Based  on  encoder 
checks  after  the  carriage  pass  these  positioning  errors  are 
consistently  less  than  +/-  .04  mm.  in  either  the  vertical  or  on-axis 
direction  (the  traverse  is  not  moved  in  the  streamwise  direction  dur¬ 
ing  a  carriage  pass) . 

Since  the  LDV  traverse  system  and  ship  model  are  attached 
separately  to  the  towing  carriage,  some  relationship  between  the 
tra'.’erse  encoders  and  the  model  coordinate  system  must  be  established 
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for  each  test.  This  can  be  done  by  placing  the  laser  beam  crossing 
point  on  model  surface  points  of  known  location.  The  encoders  can  be 
zeroed  (or  set  to  any  given  value)  at  these  points.  Encoder  readouts 
then  directly  display  distance  from  those  ship  model  points.  This 
procedure  can  be  done  best  on  black  model  surfaces.  When  low  power 
green  beams  cross  on  the  black  surface,  a  Doppler  signal  (frequency 
equals  the  effective  frequency  shift)  will  be  observed  at  the  signal 
processors  with  an  oscilloscope.  Positioning  accuracy  better  than 
+/-  .1  mm  is  possible 

These  first  three  positioning  errors  are  relatively  minor  in 
comparison  to  LDV  measurement  volume  dimensions  and  the  two  errors 
that  follow.  Combining  all  three  (they  are  independent),  the  posi¬ 
tioning  error  is  still  less  than  +/-  .11  mm. 

Major  Positioning  Errors.  The  position  encoders  actually  display 
relative  movement  between  adjacent  parts  of  the  traverse  structure. 
Bending  of  the  traverse  structure  under  its  own  weight  or 
hydrodynamic  forces  does  not  change  the  relative  position  of  adjacent 
parts.  However,  it  can  affect  the  relative  position  of  the  LDV  meas¬ 
urement  volume  and  the  independently  mounted  ship  model. 

As  traversing  shifts  the  significant  weight  of  the  optics  bread¬ 
board  and  support  members,  the  whole  traverse  structure  undergoes 
bending  changes.  Investigations  have  shown  that  most  of  the  bending 
occurs  near  the  top  of  the  structure  and  is  associated  with  the  vert¬ 
ical  traverse  mechanism  and  the  attachment  to  towing  carriage 
members.  Bending  angles,  though  small,  are  transformed  into 
significant  displacements  by  the  distance  between  the  source  of  bend¬ 
ing  and  the  point  of  interest  (the  LDV  measurement  volume). 

This  bending  is  a  fairly  repeatable  function  of  the  absolute 
traverse  location  and  therefore  positioning  corrections  can  be  made. 
These  errors  were  documented  (Figures  16  and  17)  by  careful  experi¬ 
mentation  with  a  stationary  towing  carriage.  In  certain  extreme 
cases  the  positioning  error  can  be  greater  than  10  mm.  Appendix  6 
formulates  how,  when  necessary,  these  errors  can  be  reduced  by 
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correction  factors  that  are  incorporatated  into  the  data  taking 
computer  software.  Figures  16  and  17  also  show  the  positioning  error 
reductions  possible  when  correction  factors  are  used. 

Traverse  structure  bending  also  results  from  the  two  main 
hydrodynamic  forces  exerted  on  the  submerged  strut.  They  are  drag 
and  side  forces.  Both  are  a  function  of  towing  speed,  flow  angle  of 
attack  (for  strut),  and  the  assymmetric  probe  box  fairing.  Thus, 
they  are  also  a  function  of  strut  position  in  the  flow  field  as 
modified  by  the  towed  model.  Figure  18  shows  strut  position  changes 


as  a  function  of  towing  speed.  These  measurements  were  actually  made 
on  the  strut  (as  far  down  as  possible  while  still  above  the  water 
surface),  but  they  should  closely  approximate  measurement  volume 
deflections.  Even  under  the  freestream  conditions  of  figure  18,  the 
repeatability  of  deeply  submerged  strut  deflections  is  disappointing. 
Unstable  flow  about  the  strut  may  be  the  reason.  Appendix  0  shows 
how  traversing  calibrations  under  tow  can  correct  for  both  static  and 
some  dynamic  bending  errors. 

These  last  two  positioning  errors  (due  to  bending)  depend  on 
which  corrections  are  used  and  the  measurement  point's  distance  from 
model  reference  positions  (where  there  is  no  positioning  error). 

Total  Positioning  Error.  The  five  contributors  to  the  positioning 
error  are  independent  of  each  other.  The  total  positioning  error  is 
"slow  changing"  and  "position  dependent".  From  setup  to  setup  it 
should  have  a  mean  of  zero  and  standard  deviation  of*: 

Positioning  Error  With  No  Corrections 

Pv  =  ( ( (Y-Y^) *3.0)^  +  ( ( (Z-Z^)  *5.0)^  +  0.011^)°*^  mm.  (28) 

X  o  o 

Py  =  (  ( (Y-Y^) *12.)^  +  0.011^)°*^  mm. 

=  (  (  (Z-Zy) *45. ) ^  +  0.011^)^*^  mm. 

Positioning  Error  With  Bending  Corrections 

Px  =  (  (  (Y-Y^) *0.3)^  +  (  (  (Z-Z^) *0.6)^  +  0.011^)°*^  mm.  (29) 

Py  =  (  (  (Y-Yq)  *1 . 2)^  +  O.Oll^)'^*’’  mm. 

P^  =  (  (  (Z-Z^) *4.5)^  +  0.011^)°*^  mm. 

where  X,  Y,  Z  :  vertical,  on-axis,  streamwise  position  (m.) 

X„,Y  ,Z  :  model  reference  point  location  (m.) 


dynamic  bending  errors  are  assumed  zero  or  accounted  for 
by  model  position  referencing  during  a  carriage  pass 


Nothing  has  been  said  yet  about  the  velocity  gradient  part  of 
the  positioning  velocity  error.  The  local  velocity  gradient  is 
"position  dependent”  and  "fast  changing"  but  does  not  generally  have 
a  mean  value  of  zero.  Velocity  gradients  can  only  be  guessed  or 
estimated  from  the  experimental  results  themselves.  The  highest 
gradients  occur  near  forward  ship  model  stagnation  points  and  in  thin 
model  boundary  layers.  Users  should  be  especially  careful  to  estab¬ 
lish  nearby  model  reference  points  in  these  cases. 

Velocity  Error  Due  to  Traverse  Positioning 

E^(U)  =  ((  P^*iU/ix)^  +  (  Py*>U/>y)^  +  {  P^*iU/iz)^)°-^  /  Uq  (30) 

E^(V)  =  ((  P^*W/ix)^  +  (  Py*iV/iy)^  +  (  P2*iV/dz)^)°*^  / 

E^(W)  =  ({  P^*jW/ix)^  +  (  Py*jW/iy)^  +  (  P2*iW/Jz)^)°‘^  /  Uq 

Because  of  the  difficulty  in  guessing  or  estimation  the  velocity 

gradient  at  each  measured  point,  this  error  will  be  reported  as  posi¬ 
tioning  error  (  P^ ,  P^,  P^)  and  not  a  velocity  error.  After  the 

data  has  been  taken,  the  user  must  determine  if  the  estimated  posi¬ 
tioning  error  causes  significant  velocity  errors. 

Summary  of  Finalized  Fluid  Velocity  Errors 

All  the  errors  in  this  section  are  assumed  to  be  independent  of 
each  other  and  all  previously  defined  errors.  However,  some  of  the 
suggested  correction  techniques  reduce  several  error  magnitudes  at 
once  . 

None  of  these  errors  are  "fast  changing"  and  they  will  not  be 
reduced  by  taking  the  mean  of  multiple  samples  of  the  particle 
velocity. 

Table  6  summarizes  the  error  formulations  both  with  and  without 
the  described  correction  procedures.  Uncorrected  alignment  errors 
are  probably  the  most  important  among  this  group.  Uncorrected  posi¬ 
tioning  errors  can  also  be  important  for  either  large  velocity 


TABLE  6  -  FINALIZED  WATER  VELOCITY  ERROR  FORMULAS 


Character 
Fast  Slow 

Bssssssssassaa 


Component  Velocity  Errors 
Streamwise  Vertical  On-Axis 

BSSSa3%3SSS3SSas=SSSSBSSaS33SSSSS33SS3SSSaxSSS3S3 


5*10"^ 


5*10-4 


5*10-4 


I  * 

Uncor , 


2*10-3 

0 


2*10-3 

0 


((a*U)3  +  (b*V)3  +  (c*W)3  +  (d*0g)2)0.5/ 


Uncor . 

X 

a 

= 

.  000 

a 

3 

.005 

a 

3 

.013 

b 

= 

.005 

b 

3 

.000 

b 

3 

.013 

c 

= 

.005 

c 

3 

.005 

c 

3 

.000 

d 

.000 

d 

3 

.000 

d 

3 

.000 

Corr . 

X 

a 

= 

.000 

a 

S 

.000 

a 

S 

.000 

b 

= 

.005 

b 

3 

.000 

b 

3 

.013 

c 

.005 

c 

3 

.005 

c 

3 

.000 

d 

“ 

.001 

d 

3 

,001 

d 

3 

.002 

'’x'Py.Pz* 

(  (a 

*(Y-yo))^  + 

(b*  (Z- 

•Zo 

))2  * 

0 . 011^ )  *^  •  5  nim . 

Uncor . 

X 

a 

3 

3.0 

a 

3 

12.0 

a 

3 

0.0 

b 

s 

5.0 

b 

3 

0.0 

b 

3 

45. 

Corr . 

X 

a 

3 

0.3 

a 

S 

1.2 

a 

3 

0.0 

b 

s 

0.6 

b 

3 

0.0 

b 

3 

4.5 

E  :  standard  deviation  of  velocity  error  divided  by  U^^ 

P  :  standard  deviation  of  position  error  in  units  of  ".11111." 

'Jq  towing  speed 

X,  Y,  Z  ;  vertical,  on-axis,  streamwise  position  (m.) 

^o'^o'^o  •  n’Odel  reference  point  location  (m.) 


gradients  or  large  traversing  ranges. 


VELOCITY  COMPONENT  ERROR  SUMMARY 

Tables  4,  5,  and  6  summarize  the  error  calculation  formulas 
developed  in  this  section.  In  the  case  of  a  fluctuating  velocity 
field,  generally  only: 


'  I"/-'.  <■  *■- 


V  ■.  s  . '  *•  >  *  ■- 


fc"-  ,  ' 


-L  •  tpi ....  signal  Noise  Error 

2.  E^.... Fringe  Spacing  Determination  Error 

3.  E^ . . . . Ve locity  Fluctuation  Error 

4.  E^g .. .Velocity  Bias  Error 

5.  Ep^g.. Fringe  Bias  Error 

6.  E„,_. .Filter  Bias  Error 

I?  1  D 

7.  E_ . . . . Ve loci ty  Component  Direction  Error 

a 

can  be  significant.  In  non-fluctuating  velocity  fields,  the 
calculated  bias  and  velocity  fluctuation  errors  are  zero  except  when 
FS,  and  low  pass  filter  settings  mean  a  velocity  is  totally 

unmeasureable .  However,  despite  the  calculation,  velocity  bias  can 
exist  in  a  steady  flow  field  if  a  sharp  local  velocity  gradient 
causes  the  LDV  measurement  volume  to  see  significantly  different 
velocities  at  the  same  "position". 


INITIAL  OPERATING  PARAMETER  SELECTION 

This  section  provides  a  guide  for  making  an  intelligent  estimate 
of  the  velocities  that  will  be  measured  in  an  LDV  experiment.  It 
then  helps  pick  operating  parameter  options  appropriate  to  those 
expected  velocity  component  ranges.  Finally  a  computer  program  is 
available  to  estimate  the  size  of  various  error  sources  as  a  function 
of  those  velocity  ranges  and  operating  parameters.  One  or  more 
iterations  through  this  program  can  yield  a  good  starting  point  of 
LDV  operating  parameter  choices  plus  an  idea  of  what  experimental  er¬ 
rors  to  expect. 

MEAN  VELOCITY  COMPONENT  RANGES 

The  primary  determinant  of  velocity  component  range  is  the 
experiment's  towing  speed.  This  section  will  only  discuss  component 
ranges  that  are  non-dimensionalized  by  this  speed. 

LDV  error  analysis  calculates  errors  in  each  of  the  three 

component  mean  measurements.  In  a  strict  sense,  each  component  error 

is  a  function  of  all  three  mean  velocity  components.  However,  the 

long  thin  measurement  volume  ellipsoid,  enables  very  accurate  error 

estimation  for  each  component  based  on  only  two  velocity  components. 

Some  of  these  estimates  break  down  for  velocity  directions  that  are 

less  than  15°  from  the  on-axis  direction.  Vertical  and  streamwise 

component  errors  are  primarily  functions  of  both  the  vertical  and 

streamwise  mean  velocities,  "U"  and  "V".  The  on-axis  component  error 

is  primarily  a  function  of  its  own  mean  velocity,  "W"  and  the 

2  2  0  5 

cross-measurement  volume  component  "T"  or  (U  +  V  )  . 

Figure  19  displays  velocity  "Scatter"  plots  for  two  different 
t  w 1 'Kj  tank  LDV  experiments.  There  are  two  companion  plots  from  each 
•  riment.  The  first  shows  experimental  mean  velocity  components 
:i  '■  i.'  I  for  streamwise  and  vertical  component  error  calculations.  The 
-,'.d  displays  components  for  the  on-axis  error  calculation.  Ap- 
p*’nlix  h  ,'ontains  a  more  complete  compilation  of  these  Scatter  plots 
from  LOV  experiments. 

In  estimating  mean  velocity  component  ranges,  the  LDV  user  may 
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F¥op«l«r  kiflow  Survey 
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—  With  Propetler 


Figure  19  Two  Examples  of  Scatter  Plots 


utilize  the  results  of  an  experiment  in  Appendix  H  that  has  some 
similarity  to  his  own.  Pitot  tube  results  from  past  wake  surveys  are 
another  reference  source.  Whether  either  of  these  sources  or  just 
experience  and  intuition  are  used,  the  LDV  user  needs  an  estimated 
minimum  and  maximum  mean  value  of  the  four  velocity  components:  "U", 
"V",  "W",  and  "T". 

VELOCITY  FLUCTUATION  MAGNITUDE 

Each  velocity  component  has  its  own  fluctuation  magnitude 
reflected  in  the  component  standard  deviation  at  a  measurement  point. 
In  almost  all  situations,  the  fluctuation  magnitudes  will  be  nearly 
equal  for  each  component. 

Measured  LDV  standard  deviations,  which  also  reflect  the  signal 
noise  level,  for  each  component  may  not  be  equal.  Figure  20  shows 
"Histograms"  for  the  three  measured  components  in  a  ship  bow  flow 
field  where  their  was  essentially  no  fluctuation  of  any  velocity 
component.  The  width  of  the  "Histograms"  is  due  to  each  component's 
signal  noise  level.  With  proper  operating  parameters,  the  vertical 
velocity  component  will  usually  be  the  most  noise-free  of  the  three 
measured.  Noise  standard  deviations  of  about  0.5%  of  the  tow  speed 
are  common. 

Figure  21  shows  the  distribution  of  the  vertical  velocity 
component  starrlard  deviation,  for  two  experiments.  The  bow  flow 
example  indic.tes  the  complete  absense  of  measured  turbulent  velocity 
fluctuations.  The  propeller  inflow  survey  for  a  very  clean  destroyer 
hull  shows  both  non- turbulent  potential  flow  and  turbulent  flov? 
entering  the  propeller.  Appendix  I  displays  more  examples  of  these 
bar  plots. 

Since  error  estimation  is  very  much  a  function  of  the  velocity 
fluctuation  magnitude,  the  LDV  user  should  carefully  choose  one  or 
two  values  that  he  wants  to  input  to  the  error  estimating  computer 
program.  The  bar  graphs  presented  in  Appendix  I,  other  experience, 
and  intuition  may  all  help  determine  the  choice. 


LDV  OPERATING  PARAMETER  CHOICES 


Table  7  is  a  compilation  of  information  previously  listed  in 
tables  1  and  2.  It  lists  the  LDV  operating  parameters  and  the  errors 
on  which  they  impact.  Parameters,  for  which  there  is  really  no  user 
choice  to  be  made,  are  indicated.  The  rest  are  discussed  below. 

Frequency  Shift 

Experience  has  shown  that  only  a  few  frequency  shifts  make  sense 
at  possible  model  towing  speeds.  Examples  of  these  are: 

Streamwise,  "U":  0.0,  0.2,  0.5  MHz 
Vertical,  "V":  0 . 2 ,  0 . 5 ,  1 . 0 ,  2 . 0  MHz 
On-Axis,  "W" :  0.2,  0.5,  1.0,  2.0  MHz 

Figures  22a,  22b,  and  22c  provide  more  specific  guidance  in  the 
choice  of  frequency  shift  for  each  component.  Curves  in  velocity 
space  indicate  the  velocities  for  which  the  indicated  frequency 
shifts  will  produce  8  (or  16)  fringe  crossings  for  particles  travel¬ 
ing  through  the  measurement  volume  center.  In  other  words  velocities 
for  which  measurement  is  marginally  possible.  On  one  side  of  the 
line  are  velocity  component  combinations  that  are  measureable  (22a  - 
to  the  right;  22b  -  above;  22c  -  left)  at  the  labeled  frequency 
shift  and  on  the  other  side  of  the  line  those  that  are  not  measure- 
able  . 

Depending  on  the  mean  velocity  component  ranges  and  fluctuation 
magnitude  (discussed  above) ,  the  LDV  user  should  be  able  to  narrow 
his  frequency  shift  choices  to  one  or  two  values  for  each  component. 
Note  that  unworkable  (too  low)  frequency  shifts  will  result  in  the 
calculation  of  very  high  fringe  bias  errors,  Ep^-g.  Frequency  shifts 
that  are  higher  than  necessary  have  higher  than  needed  noise  error 

Frequency  Filters 

From  velocity  component  range,  frequency  shift,  and  fringe 
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spacing,  the  LDV  user  can  calculate  the  signal  frequencies  that  the 
signal  processors  will  see.  The  formula  for  such  a  calculation  is: 

F  =  FS  +  U/Sy  (31) 

where  FS:  Frequency  shift  (MHz) 

U  :  measured  velocity  component  (m/s) 

Sy :  component  fringe  spacing  (  qm . ) 

Note  that  a  frequency  range  calculation  is  necessary  for  each 
frequency  shift  of  each  velocity  component. 

The  high  pass  filter  setting  on  the  signal  processor  is  not  too 
important  and  can  be  set  safely  below  the  frequency  range.  The  low 
pass  filter  is  important  in  eliminating  some  signal  noise,  the  LDV 
user  should  choose  one  or  possibly  two  settings  that  are  as  close  as 
possible  to  the  maximum  signal  frequency.  The  lower  the  setting,  the 
more  noise  elimination.  If  the  setting  is  too  low,  error  calcula¬ 
tions  of  filter  bias,  will  be  large  and  at  least  some  LDV  sig¬ 

nals  will  be  noticeably  attenuated. 

Input  Signal  Gain 

This  parameter  is  set  during  the  experiment  when  real  LDV  sig¬ 
nals  can  be  observed.  It  sometimes  adjusted  during  the  course  of 
measurements.  Noise  error, calculations  assume  that  it  is  adjusted 
so  that  the  noise  between  signal  bursts  has  a  peal^  to  pea]<  amplitude 
of  lass  than  50  mV.  (or  half  the  signal  processor  trigger  levels). 


Ampl i tude  Limit 

This  a<liustment  shoulh  bo  turned  off.  Experience  has  shown  that 
measured  velocities  come  only  from  seeded  particles.  If  particle  lag 
error,  E^  estimates  are  too  large  for  the  generally  used  particles 
(Figure  n)  ,  then  nst^  of  this  adjustiment,  can  be  considered.  However, 
a  surer  solution  is  to  buy  smaller  seed  particles  or  classify  the 
seed  on  hand  into  a  smaller  size  range.  This  should  normally  not  be 
nect-ssa  r  y . 

Amf.tlitude  limit  does  not  enter  into  the  error  calculations 


except  to  perhaps  change  the  assumed  mean  particle  diameter  of  6  m. 
used  in  the  calculation  of  Ej^.  A  key  drawback  in  using  the  adjust¬ 
ment  is  that  there  is  no  easily  determined  relationship  between  set¬ 
ting  and  maximum  measured  particle  size. 

Signal  Processor  Mode 

It  is  recommended  that  "TBC"  or  total  burst  count  mode  be 
employed  on  all  signal  processors.  This  allows  the  velocity  bias 
correction  (based  on  total  cycles  per  burst)  to  be  done.  This  mode 
also  permits  bad  data  elimination  by  comparing  primary  and  secondary 
fringe  crossing  timings  as  assumed  in  signal  noise  error,  calcula¬ 
tions  . 

"TBM"  or  total  burst  mode  is  the  only  other  mode  that  allows 
velocity  bias  correction  but  it  does  not  allow  timing  comparisons. 

FRINGE  COUNT, 

It  is  recommended  that  either  8  or  16  be  chosen  for  all  signal 
processors.  If  an  value  lower  than  8  is  chosen  then  the  superior 
5  to  8  timing  comparison  cannot  be  done.  The  SNR  value  used  in  sig¬ 
nal  noise  error, E^  calculations  assumes  such  a  comparison.  Also,  the 
higher  the  fringe  crossing  count,  the  greater  the  possibility  that 
random  noise  spikes  may  interrupt  proper  counter  processor  threshold 
crossings.  Such  data  would  be  invalidated  (lowering  data  rates)  or 
r>_'Corded  erroneously. 

Note  that  fringe  count  and  frequency  shift  for  each  component 
shoul'i  be  chosen  together.  Frequency  shifts  must  be  high  enough  to 
allow  ill  expected  velocities  to  cross  at  least  M^  fringes  in  the 
mei  ■•.u  r  em^' n  t  volijca.  The  best  "FS"  -  combination  trys  to: 

1.  Minimizes  the  component's  noise  error,  E^ 

2.  Avoids  large  fringe  bias  errors,  Ep^^^ 

Pui-s  rv-),'.  m.ixinuim  signal  frequency  very  close  to  a  low 
n  1  s f  1  1  tf-  r  choice. 


Timing  Comparison  Accuracy 

A  setting  of  1%  is  recommended  for  best  erroneous  data  removal. 
The  SNR  value  used  in  signal  noise  error, calculations  assumes  this 
1%  comparison  value  (Table  4). 

Velocity  fluctuation  frequencies  are  generally  much  lower  than 
necessary  to  cause  real  fringe  crossing  frequency  changes  during  a 
single  burst  timing.  If  the  user  is  unsure,  a  comparison  can  be  done 
between  (eq.#3)  and  the  expected  velocity  fluctuation 
frequency . 

Total  Number  of  Measurements,  "N" 

A  recommended  value  of  768  measurements  (256  per  component)  is 
suggested  unless  a  combination  of  tight  accuracy  requirements  and 
high  velocity  fluctuation  magnitudes  warrant  more.  Note  that  or 
E^  errors  have  to  be  major  errors  for  an  increase  in  "N"  to 
significantly  affect  the  total  error.  Also,  note  that  the  use  of  "N" 
in  error  calculations  assumes  that  each  measurement  is  independent. 
Experiment  data  rates  faster  than  the  frequency  of  velocity  fluctua¬ 
tions  do  not  supply  completely  independent  measurements. 

Five  Data  Corrections 

Previously  described  data  corrections  for  component  direction, 
flow  d  i  sturbanc.'o ,  frequency  shift,  velocity  bias,  and  positioning  er¬ 
rors  are  al.v.ys  recommend(>d .  Note  that  though  these  corrections  are 
done  during  data  analysis,  they  rely  on  certain  measurements  that 
must  be  taken  during  the  course  of  the  experiment. 

Erroneous  Data  Elimination 

The  use  of  Chauvenf>ts'  Criterion  is  recommended  in  all  locations 
excerpt  where  velocity  component  distributions  are  expected  to  be  far 
from  guissian.  Meisurement  locations  that  are  crossed  by  the 
fluctuating  edge  of  a  turbulent  boundary  layer  or  wake  are  examples 
of  when  not  to  us-'  any  data  - '  1  i  m  i  na  t  i  on .  Note  that  such  a  determina¬ 
tion  requires  cireful  review  of  the  velocity  component  histogrrims  and 
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still  IS  usually  a  judgement  call  at  best. 

ERROR  CALCULATION  COMPUTER  PROGRAM 

Two  versions  of  an  LDV  measurement  error  program  "ERRORl"  were 
written.  Both  computer  programs  are  in  "Fortran"  and  use  the  error 
calculation  formulas  of  this  report  to  calculate  a  total  error  for 
each  velocity  component.  They  differ  only  in  the  sophistication  of 
the  output  graphics  and  the  computer  on  which  they  will  run.  All 
individual  errors  are  assumed  independent  and  combined  on  that  basis. 

The  simplest  ERRORl  version  runs  on  the  POP  11-23  computer  used 
in  LDV  data  collection.  It  outputs  error  number  tables.  Error 
contour  plots  (in  velocity  space)  are  available  from  the  second  ver¬ 
sion  written  for  a  HDL-VAX  computer  generally  used  for  LDV  data 
graphics.  The  key  to  the  VAX  implementation  is  the  availability  of 
the  DTSSPLA  plotting  package  (version  9.0). 

Velocity  range  and  LDV  operating  parameters  are  entered 
interactively  by  the  ERRORl  user.  The  user  also  chooses  which  errors 
are  to  be  calculated  and  summed.  This  allows  looking  at  either 
individual  errors  or  the  total  component  error.  Before  any  errors 
are  calculated,  the  user  is  presented  with  two  tables  listing  the 
choices  he  has  made.  This  is  a  last  chance  for  changes  before  calcu¬ 
lation  proceeds.  Figure  23  displays  example  input  tables  as  they  are 
output  by  ERRORl. 

ERRORl  does  the  error  calculations,  outputs  the  results,  and 
then  allows  the  user  to  change  the  velocity  ranges,  operating  parame¬ 
ters,  errors  calculated,  or  all  three.  The  calculations  are  redone 
as  often  as  the  user  changes  his  inputs.  The  user  must  keep  track  of 
which  are  the  largest  crrcrs,  because  only  by  improving  them  can  the 
total  error  be  significantly  reduced.  In  general,  fluctuating 
velocity  fields,  only; 


,  and  E. 


can  be  significant.  Of  these  E^  is  fixed  and  E  only  responds  to  the 

I  cl 

directional  dat.3  correction.  Tn  non- f  1  uc:  t  ua  t  i  ng  velocity  fields,  the 


U  Comp. 

V  Comp. 

W  Comp. 

T  Vel 

Minimum  (m/s) 

0.30 

-1.00 

-1.00 

0.30 

Maximum  (m/s) 

3.00 

1 .00 

1 .00 

3.10 

Std  Dev  (/Uo) 

0.07 

0.07 

0.07 

0.07 

Freq  Shift  (MHz) 

0.000 

0.500 

0.500 

LP  Filter  (MHz) 

2.000 

1 .000 

1 .000 

Fringe  Count  Mf 

8 

8 

16 

Meas.  Count  N 

256 

256 

256 

Comb.  Correction 

X 

X 

X 

V  Bias  Correction 

X 

X 

X 

Tow  Speed  =  3-000  (m/s) 


Signal  Noise  Error . E<n>  X 

Time  Digitization  Error . E<d>  X 

Fringe  Spacing  Error . E<f>  X 

Frequency  Shift  Error. . E<s>  X 

Particle  Lag  Error . E<L>  X 

Velocity  Fluctuation  Error. ..E<v>  X 

Filter  Bias  Error . £<F1B>  X 

Carriage  Speed  Error . E<c>  X 

Flow  Disturbance  Error . E<1>  X 

Velocity  Bias  Error . E<VB>  X 

Fringe  Bias  Error.... . E<FrB>  X 

Vel.  Comp.  Direction  Error. ..E<a>  X 


XXXXXX/XXXXXX/XXXX.XXXXX.XXX 

xxxxxx/xxxxxx/----.  --- 

xxxxxx/xxxxxx/----.-----.--- 


Figure  21  Errorl" Input  Tables 


calculatfd  bias  and  velocity  fluctuation  errors  are  zero  except  when 
FS,  and  low  pass  filter  settings  mean  that  a  velocity  is  totally 
unmeasureable .  However,  despite  the  calculation,  velocity  bias  is 
possible  if  a  sharp  local  ve?locity  gradient  causes  the  LDV  measure¬ 
ment  volume  to  see  significantly  different  velocities  at  the  same 
"position".  Note  positioning  error  is  not  addressed  by  ERRORl. 

Example  outputs  for  the  throe  velocity  component  errors  appear 
in  Figure  24a  and  24b.  The  simplier  output  (POP  11-23  version)  and 
the  contour  output  { VAX -D I SSPLA  Version)  are  displayed  side  to  side 
for  the  same  inputs.  The  simpler  output  ve^rsion  is  analogous  to  the 
contour  plot.  Rut  instead  of  contours,  error  values  on  a  regular 


grid  (in  velocity  space)  appear.  Note  that  the  plot  axes  are  exactly 
those  used  previously  for  Scatter  plots  (Figure  19)  and  for  Velocity 
measureabi 1 i ty  (Figure  22).  In  fact  a  superposition  of  the  error 
contour  plot  and  a  corresponding  Scatter  plot  is  suggested  to  visual¬ 
ize  and  keep  track  of  experimental  errors. 

The  user  should  carry  away  from  his  session  with  ERRORl,  a  set 
of  error  contour  plots  corresponding  to  his  experiment's  optimum 
(smallest  total  error)  operating  parameters.  Actually  there  may  be 
more  than  one  such  plot  per  component  because  of  uncertainties  in 
mean  velocity  ranges  and  velocity  fluctuation  magnitudes. 

Several  notes  about  the  ERRORl's  characteristics  are  helpful  in 
making  efficient  use  of  the  program.  The  first  is  that  though  the 
program  is  conveniently  available  on  the  small  POP  11/23,  it  just 
barely  fits  into  memory  and  a  special  linking  procedure  is  required 
to  produce  a  executable  file.  In  addition  whenever  velocity  or 
fringe  bias  errors  are  calculated,  the  error  calculation  procedure 
takes  3  minutes  or  more  to  complete.  It  takes  a  few  seconds  if  these 
two  errors  are  not  calculated.  The  VAX  version  (when  the  system  is 
not  overloaded  with  users)  is  several  times  faster  in  the  calcula¬ 
tions,  but  each  contour  plot  takes  about  a  minute  to  plot  on  the 
terminal  screen.  As  a  printed  message  explains,  unchanged  inputs  do 
not  have  to  be  continually  reentered  for  each  error  calculation.  Fi¬ 
nally  both  program  versions  keep  track  of  what  inputs  have  changed 
from  one  error  calculation  to  the  next.  If  nothing  has  been  changed 
to  affect  a  given  error,  the  program  simply  uses  the  previously 
calculated  value. 
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MONITORING  EXPERIMENTAL  ERRORS 


During  LDV  experiments,  a  careful  user  should  monitor  to  some 
degree  the  errors  of  the  measurements  he  is  taking.  Also,  he  should 
complete  all  calibrations  and  error  checks  that  will  be  needed  later 
during  data  analysis.  This  section  outlines  the  procedures  involved. 

PARAMETER  SETTING  AND  ERROR  TRACKING  DURING  EXPERIMENT 

Position  correction  schemes  are  described  in  Appendix  G  based  on 
either  static  (in  drydock)  or  dynamic  (towed  LDV  system)  measure¬ 
ments.  Dynamic  corrections  increase  in  importance  for  faster  and 
deeper  LDV  struts.  In  addition  to  position  corrections,  the  "in-tow" 
condition  provides  the  best  time  to  establish  the  traverse  encoder 
location  of  some  model  reference  point(s).  Both  LDV  traverse  and 
ship  model  strut  bending  can  change  with  tow  speed.  Either  way  the 
establishment  of  model  reference  points  is  affected. 

The  user  should  take  care  in  applying  the  position  correction 
scheme.  As  a  quick  check  of  proper  usage; 

1.  Traverse  the  full  "Y'"  range  (  "I.  m)  using  the  LDV 

software  and  corrected  positions;  compare  the  encoder 

readout  change  to  the  corrected  position  change;  the 

encoder  readout  should  be  on  the  order  of  10  mm  larger. 

2.  Traverse  the  full  "Z"  range  (  ~0.5  m)  using  the  LDV 

sof^ware  and  corrected  positions;  compare  the  encoder 

readout  change  to  the  corrected  position  change;  the 

encoder  readout  should  be  on  the  order  of  25  mm  larger. 

3.  Compare  software  corrected  "X"  (vertical)  positions  to 
encoder  readouts:  the  encoder  readout  should  always  be 
greater  than  or  equal  to  the  corrected  position  (i.e., 
higher)  but  by  no  more  than  1  or  2  mm. 

Each  experiment  should  include  some  initial  carriage  runs  at  tow 
speeds  that  will  be  used  during  measurements.  By  observing  the 
doppler  bursts  on  an  oscilloscope  and  the  computer  collected  data, 
the  validity  of  velocity  range  and  operating  parameter  (arrived  at  in 
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the  previous  section)  can  be  checked.  The  user  should  make  sure; 

1.  Enough,  but  not  too  many,  fringe  crossings  per  Doppler 
burst 

2.  The  next  highest  low  pass  filter  setting  does  not 

noticeably  increase  burst  amplitudes  or  data  rates 

3.  Doppler  bursts  are  pretty  clean  with  nearly  sinusoidal 
variation  with  fringe  crossings 

4.  Measured  mean  velocity  components  are  not  out  of  the  as¬ 
sumed  range 

5.  Measured  vertical  component  standard  deviation  is 

consistent  with  the  assumed  velocity  fluctuation  magnitude 

6.  The  relationship  between  data  rate  and  dispersed  seed 

particles  is  consistent  with  past  experience 

Past  experience  has  shown  on  carriage  2  that  an  average  of  1  table¬ 
spoon  of  seed  per  pass  can  yield  data  rates  for  all  components  of 
>  15*Uq  (sec  ^  where  is  in  m/s).  The  streamwise  and  on-axis 

component  data  rates  are  significantly  slower  than  the  attainable 
vertical  component  data  rate. 

If  the  previous  "ERRORl"  calculations  are  not  applicable  to  the 
final  parameters  or  velocity  ranges  found,  then  new  calculations 
should  be  done.  With  these  error  contours,  the  user  can  start  real 
measurements.  He  can  locate  by  hand  some  of  his  data  on  the  error 
contour  plots.  In  doing  so,  he  should  be  on  the  look  out  for: 

1.  Mean  velocities  exceeding  plot  ranges 

2.  Data  points.  occurring  in  regions  of  unacceptably  high 
error 

3.  Vertical  Velocity  standard  deviations  exceeding  the  as¬ 
sumed  velocity  component  fluctuation  magnitude 

The  occurrenc^^  of  any  of  these  may  require  operating  parameter 
changes  to  minimize  measurement  error.  Often  the  need  for  parameter 
changes  is  dramatically  indicated  by  component  data  rates  falling  to 
zero  or  measured  standard  deviations  becoming  unexpectedly  large. 
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RECORDS  FOR  FUTURE  DATA  ANALYSIS: 

Fringe  spacing  for  the  streamwise  and  vertical  components  are 
determined  by  making  measurements  on  a  spinning  wheel  surface  that 
travels  at  a  very  precise  speed.  An  error  estimate  for  this  determi¬ 
nation  was  made  previously  in  the  discussion  of  E^ .  However,  certain 
assumptions  in  this  estimate  may  not  be  justified  (a  perfectly 
concentric  rotating  wheel,  equal  fringe  spacing  throughout  the  meas¬ 
urement  volume,  etc.).  So,  instead  of  one  fringe  spacing  determina¬ 
tion  (always  needed  to  begin  taking  data),  the  user  should  take  four 
or  more.  Preferrably  the  determinations  should  be  independent  setups 
of  the  calibration  wheel  and,  if  possible,  at  both  the  beginning  and 
end  of  any  optical  alignment. 

The  description  of  the  velocity  component  directional  error 
included  a  description  of  how  LDV  measurements  taken  with  no  ship 
model  present  could  by  used  to  wholly  or  partially  correct  for 
several  errors.  This  procedure  which  consists  of  two  or  three  car¬ 
riage  passes  worth  of  data  should  be  repeated  for  every  optical 
alignment  during  the  experiment.  Table  8  displays  data  obtained  from 
such  a  procedure.  The  non-zero  values  for  the  on-axis  and  vertical 
components  are  assumed  due  to  their  misalignment  with  the  streamwise 
direction.  The  difference  of  the  streamwise  velcity  from  the  tow 
speed  is  assumed  tro  be  due  to  fringe  spacing  measurement  error,  flow 
disturbance  by  the  LDV  strut,  or  both. 

Finally  a  data  repeatability  test  should  be  run  at  some  time 
during  the  experim.ent.  This  means  that  data  should  be  repeated  at  a 
few  locations  (at  least  20  for  every  optical  alignment)  to  allow 
la^-.or  evaluation  of  raeasr,?  ement  repeatability,  as  in  Figure  25. 
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FINAL  ERROR  ESTIMATION 


A  recommended  three  component  LDV  data  analysis  procedure  is 
detailed  in  Reference  13.  This  procedure  contains  many  steps  to  cor¬ 
rect  for  or  evaluate  measurement  errors.  This  section  will  repeat 
the  description  of  those  data  analysis  steps  that  concern  measurement 
errors . 

HAND  CALCULATIONS  OF  CALIBRATION  DATA 

The  previous  section  on  experimental  procedures  explained  how 
both  freestream  LDV  measurements  and  fringe  spacing  calibrations  were 
to  be  taken.  Assuming  this  information  was  properly  recorded,  it  is 
relatively  simple  to  extract  the  needed  calibration  factors. 

Data  like  that  displayed  in  table  8  can  be  averaged  to  calculate 
the  correction  factors  (Ay,  A^,  A^,  By,  B^,  and  B^)  as  described  in 
the  Velocity  Component  Directional  Errors  section.  As  previously 
mentioned,  the  "B"  factors  are  of  marginal  importance  and  can  usually 
be  assumed  to  be  zero. 

Freestream  LDV  measurements  of  velocity  component  standard 
deviations  are  unique  for  each  experiment's  optical  alignment  and 
operating  parameters.  The  noise  error  for  each  component,  E^,  is  ac¬ 
curately  estimated  by  the  averaging  component  standard  deviations  ob¬ 
tained  in  the  steady  freestream  flow. 

The  four  or  more  fringe  spacing  calibrations  done  for  each  opti¬ 
cal  alignment,  should  be  combined  to  estimate  a  mean  and  standard 
deviation  of  this  measurement.  The  reliability  of  these  estimates  is 
directly  related  to  the  number  of  independent  calibrations  done  dur¬ 
ing  the  experiment.  A  separate  calculation  is  required  for  stream- 
wise  and  vertical  calibrations.  Since  the  on-axis  fringe  spacing  is 
calculated  from  the  streamwise  spacing,  the  means  and  standard  devia¬ 
tions  of  the  two  are  related  (eq.  #7  &  8 )  . 

COMPUTER  PROGRAMS  USEFUL  IN  DATA  ERROR  ANALYSIS 

This  section  describes  the  computer  programs  used  in  the  data 
analysis  procedure^^  that  deal  with  measurement  errors.  Listings  of 
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the  programs  discussed  appear  in  References  14  along  with  other  data 
analysis  and  data  display  programs.  Here,  each  program  is  discussed 
separately  in  terms  of  its  purpose,  user  interaction,  and  program 
output . 

NORMLZ ;  Data  NORMaLiZation  Program 

NORMLZ  takes  any  TSI  Software  Format  data  table  (e.g.  see  Table 
8)  and  performs  certain  user  controlled  corrections  and/or  normaliza¬ 
tions.  This  includes  the  component  directional  error  correction  and 
adjustment  of  position  coordinates  by  an  additive  constant  and 
non-dimensional izing  length  scale.  The  program  is  meant  to  run  on 
the  PDP  11-23  data  taking  computer  or  equivalent. 

User  Interaction.  The  user  must  create  an  input  file:  NORMLZ.DAT  to 

contain  the  direction  correction  factors  (A,,,  A„,  A„,  ,  B...  B,,)  and 

U  V  W  U  V  w 

position  coordinate  additive  constants  and  non-dimensional izing 
lengths.  The  exact  format  of  this  file  is  shown  in  reference  13  and 
14.  While  running  the  program,  the  user  is  asked  to  verify  if  each 
of  these  adjustments  are  actually  to  be  used.  If  the  user  enters  no, 
then  that  particular  correction  or  adjustment  is  skipped. 

The  program  is  written  to  deal  with  a  limited  number  of  data 
adjustments.  It  can  easily  be  rewritten  to  cover  a  variety  of  situa¬ 
tions  that  can  occur  during  data  aquisition  (e.g.  erroneous  fringe 
spacings  or  recorded  carriage  speeds). 

Program  Output.  The  program  outputs  a  data  file  with  table  titles 
and  data  formats  exactly  the  same  as  the  input  data  file.  Only  the 
position  and/or  velocity  values  have  been  adjusted.  A  few  table  ti¬ 
tles  are  also  changed  to  reflect  position  coordinate 
non-d imens i onal i za t i on  if  it  has  been  chosen. 


SCATRl  &  SCATR2:  Mean  Velocity  Component  SCATteR  Plots 

13 

These  two  programs  take  "Plot-able"  LDV  data  files  and  mark 
all  data  points  on  a  graph  whose  axes  are  two  velocity  components. 
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SCATRl  plots  data  for  "U"  (streamwise  component)  and  "V"  (vertical 
component)  axes.  SCATR2  plots  data  for  ”T"  (  [U^+V^]*^  )  and  "W" 

(on-axis  component) .  The  program  is  written  in  Fortran,  but  includes 
calls  to  the  plotting  package,  DISSPLA.  The  program  has  been  run  on 
the  HDL  VAX  computer. 

This  program  is  useful  in  conjunction  with  error  contour  plots 
output  by  ERRORl  (described  in  Error  Calculation  Computer  Program 
section) .  Together  they  show  an  estimate  of  the  range  of  measurement 
errors  for  any  component.  Also,  the  plots  may  be  used  in  establish¬ 
ing  velocity  component  ranges  for  future,  similar  experiments.  The 
required  frequency  shifts  and  other  LDV  operating  parameters  can  be 
established  when  used  with  Figures  22a,  22b,  and  22c. 

User  Interaction.  The  user  copies  the  data  file  of  interest  to  a 

file  named  "XWAKE.DAT"  and  then  just  runs  the  program  on  a  Tektronix 

or  other  graphics  display  terminal.  The  axes  are  labeled  in  velocity 

components  divided  by  the  towing  speed,  U^.  If  the  user  wishes  other 

than  the  default  axes'  bounds,  than  he  need  only  change  the  "Call 

15 

GRAF (....)"  program  statement. 

Program  Output.  The  program  outputs  a  graph  with  "X"s  marking  the 
location  of  each  data  point  in  the  LDV  data  file  (Figure  19  and  Ap¬ 
pendix  H) .  The  axes  are  automatically  set  and  labeled.  An  error 
should  result  if  data  points  are  far  outside  the  plot  boundaries. 

BGRAFl;  Velocity  Fluctuation  Magnitude  Bar  GRAph 

This  program  looks  at  the  standard  deviations  of  any  of  the 
measured  velocity  components  as  found  on  a  standard  "Plot-able"  data 
file.  A  bar  graph  is  produced  that  shows  the  distribution  of 
standard  deviation  values  for  any  one  velocity  component.  The 
program  is  written  in  Fortran,  but  includes  calls  to  the  plotting 
package,  DISSPLA.  The  program  has  been  run  on  the  HDL  VAX  computer. 

The  plots  (especially  the  vertical  standard  deviation)  may  be 
useful  in  establishing  velocity  component  fluctuation  magnitudes  for 
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future,  similar  experiments.  Derived  values  input  to  ERRORl  when 
future  experimental  LDV  operating  parameters  are  being  determined, 
would  increase  the  reliability  of  error  calculation  results. 

User  Interaction.  The  user  copies  the  data  file  of  interest  to  a 
file  named  "XWAKE.DAT"  and  then  runs  the  program  on  a  Tektronix  or 
other  graphics  display  terminal.  The  program  asks  for  the  component 
standard  deviation  to  be  plotted.  The  upper  bounds  of  both  the 
standard  deviation  and  the  percent  occurrence  axes  are  entered  by  the 
user  after  he  is  prompted  by  the  computer  with  the  maximum  values 
that  occur  in  the  data  set. 

Program  Output.  The  program  outputs  a  bar  graph  like  that  shown 
earlier  in  Figure  20  and  in  Appendix  I. 

BGRAF2;  Data  Repeatability  Bar  GRAph 

This  program  looks  at  the  repeatability  of  the  measured  mean 
velocity  components  as  found  on  any  TSI  Software  Format  data  file. 
The  file  must  contain  some  locations  where  two  or  more  valid  measure¬ 
ments  were  made.  At  a  repeated  location,  the  difference  of  each 
component  measurement  from  the  mean  measured  value  is  computed.  A 
bar  graph  is  produced  that  shows  the  distribution  of  the  absolute 
value  of  these  differences  for  a  component  chosen  by  the  user.  The 
program  is  written  in  Fortran,  but  includes  calls  to  the  plotting 
package,  DISSPLA.  The  program  has  been  run  on  the  HDL  VAX  computer. 

The  plots  are  useful  in  demonstrating  the  quality  of  the 
velocity  measurements  and  the  steadiness  of  model  conditions  over 
different  carriage  passes,  different  days,  or  different  experiments. 

User  Interaction.  The  user  copies  the  TSI  Software  Format  data  file 
of  interest  to  a  file  named  "REPEAT.DAT"  and  then  runs  the  program  on 
a  Tektronix  or  other  graphics  display  terminal.  The  program  asks  for 
the  particular  velocity  component  to  be  plotted.  The  upper  bounds  of 
the  both  the  repeatability  and  the  percent  occurrence  axes  are 


entered  by  the  user  after  he  is  prompted  by  the  computer  with  the 
maximum  values  that  occur  in  the  data  set. 

Program  Output.  The  program  outputs  a  bar  graph  like  that  shown 
earlier  in  Figure  25. 

BGRAF3:  Streamwise  Fringe  Crossings  Bar  GRAph 

This  program  looks  at  the  frequency  of  occurrence  of  different 
numbers  of  fringe  crossings  for  the  streamwise  LDV  velocity 
component.  The  data  is  taken  from  calibration  runs  (no  ship  model) 
when  the  velocity  direction  is  essentially  streamwise.  The  plot  is 
in  the  form  of  a  histogram  or  bar  graph.  The  program  is  written  in 
Fortran,  but  includes  calls  to  the  plotting  package,  DISSPLA.  The 
program  has  been  run  on  the  HDL  VAX  computer. 

The  input  file  called  "FRINGE.DAT”  must  be  obtained  by  examining 
the  Raw  Data  files  for  a  calibration  run.  The  data  must  have  been 
recorded  with  the  streamwise  signal  processor  in  TBC  or  TBM  mode. 

User  Interaction.  The  user  cannot  use  any  of  the  normal  data  table 

files,  but  must  somehow  decode  the  octal  words  of  raw  data  files.  He 

must  extract  from  only  streamwise  words  the  total  fringe  crossing 
3  4 

number  '  .  A  Fortran  program  FRINGE  was  written  for  the  PDP  11-23 

data  taking  computer  to  do  this  job.  There  is  only  a  hard  copy 

C’tput  from  which  the  user  must  generate  a  FRINGE.DAT  file  formatted 

14 

for  input  to  BGRAF3  .  Afterwards  the  user  runs  the  program  on  a 
Tektronix  or  other  graphics  display  terminal.  The  program  asks  for 
the  particular  component  to  be  plotted.  The  upper  bounds  of  the  both 
the  fringe  crossings  and  the  percent  occurrence  axes  are  entered  by 
the  user  after  he  is  prompted  by  the  computer  with  the  maximum  values 
that  occur  in  the  data  set. 


ERR0R2:  Final  ERROR  Estimation 


.  1"^ 

This  program  takes  "Plot-able"  LDV  data  files  and  estimates 

the  component  errors  for  each  data  point.  The  calculations  are  based 

on  the  error  formulations  of  this  report.  Unlike  ERRORl,  velocity 

ranges  and  fluctuation  magnitudes  are  not  entered  by  the  user. 

ERROR2  obtains  the  information  from  the  measurements  at  each  data 

point.  The  program  is  written  in  Fortran,  but  includes  calls  to  the 

plotting  package,  DISSPLA.  The  p  ogram  has  been  run  on  the  HDL  VAX 

computer . 

This  program  gives  the  best  possible  estimate  of  the  errors 
present  in  each  velocity  component  at  each  measurement  point. 

User  Interaction.  The  user  copies  the  data  file  of  interest  to  a 

1  file  named  "LDA.DAT"  and  then  runs  the  program  on  a  Tektronix  or 

other  graphics  display  terminal.  The  user  is  asked  to  enter  LDV 
system  operating  parameters  and  the  errors  to  be  summed  (exactly  as 
in  ERRORl) .  The  user  is  also  asked  for  estimates  of  freestream 
component  standard  deviation  magnitudes  (assumed  equal  to  E^^)  .  Also, 
a  standard  deviation  of  the  repeated  fringe  spacing  calibrations  is 
entered  for  each  velocity  component  (assumed  equal  to  E^) .  The 
origin  of  both  of  these  last  two  inputs  was  discussed  earlier  in  this 
section  under  Hand  Calculations  of  Calibration  Data. 

Program  Output.  Output  is  in  two  possible  forms.  A  bar  graph  can  be 
displayed  that  plots  any  component's  total  error  magnitude  versus 
frequency  of  occurrence.  Also,  a  revised  data  table  (ErLDA.DAT)  can 
be  output  with  measured  velocity  standard  deviations  replaced  by 
estimates  of  the  total  error.  Total  error  means  the  standard  devia¬ 
tion  expected  for  measured  component  means  about  the  true  velocity 
component  mean  (as  in  ERRORl). 

The  data  table  ErLDA.DAT  can  be  quite  helpftil  in  any  LDV  data 
plot.  Reference  to  this  table  gives  error  "bars"  or  limits  for  each 
velocity  component  at  each  measurement  point. 
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REPORT  SUMMARY 

There  is  essentially  only  one  way  to  physically  set  up  the  tow¬ 
ing  tank  LDV  system.  However,  once  it  is  in  place,  the  user  is  con¬ 
fronted  with  many  operating  parameter  choices  before  taking  any  data. 

This  report  attempted  to  guide  the  user  in  making  operating 
parameter  choices  as  well  as  enable  him  to  quantify  the  errors 
present  in  his  measurements  at  various  times.  Capabilities  and  er¬ 
rors  are  first  of  all  a  function  of  the  LDV  equipment  design  and  flow 
field  characteristics.  However,  they  are  also  highly  dependent  on 
operating  parameters  chosen  by  the  user. 

The  sections  of  this  report  were  of  two  types: 

1.  Background  Information 

2.  Reference  Information 

Two  sections  initially  supplied  background  information.  The  first 
listed  the  operating  parameter  choices  for  the  Towing  Tank  system  and 
described  how  they  affect  LDV  measurements.  The  second  listed  all 
possible  measurement  errors  and  established  their  functional  rela¬ 
tionship  to  operating  parameters,  LDV  equipment  design,  and  flow 
field  characteristics.  Various  appendixes  supported  and  amplified 
the  information  given. 

Three  reference  sections  followed.  The  first  provided  guidance 
in  making  operating  parameter  choices  and  measurement  error  estima¬ 
tions  prior  to  an  experiment.  The  second  described  how  during  an 
experiment,  the  user  can  determine  if  proper  parameters  choices  have 
been  made.  The  last  section  detailed  how  final  error  estimates  can 
be  made  after  an  experiment  for  inclusion  in  written  reports. 

The  main  purpose  of  this  report  is  to  provide  experimenters  with 
the  three  reference  sections  on  how  to  prepare  for,  execute,  and 
analyze  LDV  data  with  the  minimum  of  measurement  error.  Procedures 
and  computet  programs  to  facilitate  proper  use  of  the  towing  tank  LDV 
system  are  detailed.  The  detail  of  the  background  section  and  most 
of  the  appendixes  are  for  the  user  who  wants  to  know  more  about  the 
operation  and  errors  of  LDV  systems.  In  making  unusual  measurements 
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by  modifying  the  more  or  less  standard  mean  velocity  measurement 
procedures,  this  knowledge  may  be  very  important. 
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APPENDIX  A  TRAHSMITTIHG  OPTICS  COMPONENTS 
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9102-12 

0 

0 

0 

0 

3eam  Splitting  Cube 

y2l6 

o 

0 

0 

0 

Polarization  Rotator 

9102-12 

0 

0 

0 

‘,0  mm  Beam  Splitter 

9115-1 

0 

0 

"O' 

k 


^  !**•  n‘. 


»*' -  V 


^v‘'.'.'!v^.';.  -V-.-'v'  ^ 
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APPENDIX  A  TRANSMITTING  OPTICS  COMPONENTS 


Component 

Description 


I  TSI  Model  )  Affected  !  Affected 
I  No.  I  Blue  Beams  I  Green  Beams 

B  I  C 


Beam 

Pattern 


40  MHz  Bragg  Cell 

9182-11 

Beam  Steering  Wedges 

9175 

40  MHz  Bragg  Cell 

9182-12 

Beam  Steering  Wedges 

9175 

60  MHz  Bragg  Cell 

9282-12 

Beam  Steering  Wedges 

9175 

Beam  Spacing  Adjuster 

9220 

Displacement  Prism 

9174 

Beam  Blockers 

9181-3 

22  mm  Beam  Spacers 

9113-22 

22  mm  Beam  Spacers 

91 13-22 

2.2x  Beam  Expander 

9188 

140  mm  Beam  Spacers 

9210 

Mirror  (150  mm  Dlam.) 

Special 

Mirror  (150  x  80  mm) 

Special 

Lens  (F.L.  s  600  mm) 

Special 

APPENDIX  B;  LASER  BEAM  WAVELENGTHS 


Laser  beam  wavelengths  from  transmission  in  air  are  well  known 
and  highly  repeatable.  The  two  colors  used  in  the  towing  tank  LDV 
are  specified  to  the  tenth  of  a  nanometer. 

Xgreen  514.5  nm  (in  air)  (B-1) 

^blue  “  488.0  nm  (in  air)  (B-2) 

The  light  wavelength  in  water  is  modified  by  water's  index  of  refrac¬ 
tion,  n.  The  index  of  refraction  of  water  is  most  affected  by  water 
temperature.  Data  from  reference  13  shows  this  variation  for  pure 
water  and  589.3  nm.  light. 

^water  ~  ^air^  (B“3) 

Hfjiii  "n" 

14°C  1.33348 
15°C  1.33341 
16°C  1.33333 
18°C  1.33317 
20°C  1.33299 
22°C  1.33281 
24°C  1.33262 
26°C  1.33241 
28'’C  1.33219 

If  the  uncertainty  in  the  basin  water  temperature  is  2°C  then 
the  uncertainty  in  the  light  wavelengths  in  water  is  .02%  (if  towing 
basin  water  refractive  index  follows  the  table  above) . 
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APPENDIX  C;  PARTICLE  LAG  ERROR  CALCULATIONS 

Given  that  Stokes  Law  (a  "creeping  flow"  fluid-particle  interac¬ 
tion)  applies  to  the  particle  dynamics,  the  spherical  particle  equa¬ 
tion  of  motion  is:^ 


SG  *  Vol.*  dV^/dt  =  3*T1*y'*D  *  (U-V  ) 


where  Vol.:  volume  of  particle  =  0.1667*n*D^ 
V  :  fluid  kinematic  viscosity 
U  : 


(C-1) 


fluid  velocity  magnitude 
particle  diameter 
particle  velocity  magnitude 


Sinusoidal  Velocity  Variation 

Assuming  a  fluid  velocity  with  mean  magnitude  U  and  a 
sinusoidally  varying  component  (amplitude  u^;  frequency  f  the 
particle  velocity  response  can  be  derived: 


Vp(t) 


where  Tp: 


0  : 


U  -t-  (Uq*COS  (2*TI*F^*t 


SG  *D  (18*V) 

2*n*F  *r 
u  p 


0)/(l  -t-  (2*Tt*Tp*F^j)^)°-^ 


(C-2) 


The  error  magnitude  of  interest  is  the  standard  deviation  of  the 
difference  between  Vp  and  U.  This  is  approximately: 


=  (l/TT)  *  (2*'rT*Tp*F^^)  ^  *  (uy  Uq) 


(C-3) 


Turbulent  Flow  Field  Accelerations 

Using  equation  C-1,  the  particle  lag  velocity  can  be  written: 


(n-Vp)  =  (dVp/dt)*SGp*Dp  /  (18*y) 


(C-4  ) 


102 


i) 


In  turbulent  wake  flows  estimates  of  the  acceleration,  dV  /dt,  can  be 

P 

made  for  fluctuations  of  various  length  scales.  First  we  assume  that 


dV  /dt  and  dU/dt  are  approximately  equal  (which  will  be  true  when  lag 


errors  are  small)  and  then  go  on  to  estimate  dU/dt. 

For  large  scale  turbulent  eddies  the  length  scale  is  ,  the  wake 
boundary  layer  thickness.  The  velocity  scale  is  U^,  the  towing 
speed.  Thus  the  acceleration  scale  is: 


dVp/dt  «  dU/dt  «  Uq  /S 


(C-5) 


The  smallest  scale  turbulent  eddy  that  can  be  detected  by  the 

LDV  is  one  whose  length  scale  is  equal  to  twice  the  LDV  measurement 

14 

volume  size,  Turbulent  scale  arguments  assert  that  turbulent 

velocity,  acceleration,  and  time  scales  are  a  function  of  only  the 


eddy  length  scale  of  interest  and  the  energy  production  rate 


(E  w  boundary  layer.  Thus  the  "small”  scale  accelera¬ 


tion  magnitude  are  on  the  order  of; 


dV  /dt  »  dU/dt  «  (2*D„„) 


mv' 


(C-6) 


Because  the  LDV  measurement  size  will  usually  be  smaller  than  ship 
model  boundary  layers,  the  "small  scale"  acceleration  scale  will  be 
the  larger  of  the  two  acceleration  magnitudes.  Indeed  if  it  is  not, 
we  have  no  hope  of  detecting  the  boundary  layer  at  all.  The  largest 
turbulence  induced  particle  lag  incorporates  eq.#C-6  into  its  lag  er¬ 
ror  estimation: 


=  (Uq  *  SGp  *  Dp^)  /  (18*V* 


mv' 


(C-7) 
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APPENDIX  D:  BIAS  ERRORS  FOR  FREQUENCY  SHIFTED  LDV 

The  equations  in  this  appendix  are  a  minor  modification  of  equa¬ 
tions  found  in  reference  11.  The  modifications  allow  evaluation  of 
biases  for  frequency  shifted  LDV  systems  with  "moving  fringes". 
Symbols  correspond  to  those  used  in  reference  11. 

Stationary  fringes  require  for  a  measurement  that; 

/it  >  M^*s  /  U^  (D-1) 

where  /it:  transit  time  of  particle  across  measurement  volume 
M^:  fringe  crossings  needed  for  valid  measurements 

s  :  Fringe  spacing 

:  velocity  component  perpendicular  to  fringe  planes 

Moving  fringes  (with  frequency  shift,  f^)  require  for  a  measure¬ 
ment  that: 
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shifted  LDV  systems  this  equation  becomes: 


B(U;  €') 
2 

where  W  ; 

0  : 


|U^1*W*(1  -  if  £'^*W^  <  1 

0  if  Z  1 

1  +  ((Uy*sin  0/2)^  +  /  ((U^*cos  0/2) 

beam  crossing  angle 


(D-6) 


This  equation  combines  velocity  and  fringe  bias  into  one  equa¬ 
tion.  In  an  LDV  setup  (given  *  value)  with  uniform  particle  concen¬ 
tration,  the  data  rate  for  any  particular  velocity  U  is  proportional 
to  "B" .  Thus  1/B  can  be  used  as  a  weighting  factor  to  determine  mean 
values  of  the  velocity  component  from  data  acquired  from  all 

particles  passing  through  the  measurement  volume.  These  bias  effects 
can  be  simply  separated: 

B(U;  €')  =  G^(U)  *  G£(U)  (D-7) 

where  G  :  |U  |*W 

^  ^  2  2 
G^:  (1  -  €'  *W^) 


On-Axis  Component  Bias 

The  equations  in  reference  11  and  this  appendix  can  be  directly 
applied  to  biases  in  the  streamwise  and  vertical  components  of  the 
towing  tank  LDV  system.  Symbols  subscripted:  "y"  are  in  the  on-axis 
direction,  "x"  are  in  the  direction  of  the  component  being  consid¬ 
ered,  "z"  in  the  remaining  direction.  For  biases  in  the  on-axis 
direction  however,  the  following  definitions  must  be  used: 


=  {U +  U  ^  +  (U  *sin  0/2)^)/U  ^  (D-8) 

X  z  y  y 

B(u;  €')  =  |Uy|*W*(l  -  €'^*w^)  (D-9) 

y  ;  on-axis  direction 
x,z:  streamwise  and  vertical  directions 


where 


APPENDIX  E:  VARIATION  OF  MEASUREMENT  VOLUME  CROSSING  TIME 

The  difference  between  the  two  schemes  of  velocity  bias  correc¬ 
tion  is  the  choice  of  weighting  factors.  One  scheme  uses  a  computed 
rate  of  arrival  of  particles  (proportional  to  1/T_) .  The  other 

scheme  uses  individual  measurements  of  crossing  times,  T  . .  The 

B 1 

individual  measurements  vary  around  T  according  to  the  degree  which: 

B 

1.  a  particle  trajectory  misses  the  measurement  volume 
cente  r 

2.  noise  causes  the  signal  processor  to  perceive  a 
premature  end  of  the  LDV  signal  burst 

This  variation  induces  an  additional  source  of  variation  of  weighted 
velocities  about  their  true  mean. 

An  estimate  of  the  standard  deviation  T  .  can  be  calculated  by 
assuming  that  every  particle  that  passes  through  the  measurement 
volume  ellipsoid  is  measurable  and  there  is  no  "on-axis"  velocity. 
This  is  a  conservatively  high  estimate  of  the  T^,  standard  deviation 
because  it  includes  the  very  short  Tg^  values  for  trajectories  that 
barely  penetrate  the  ellipsoid  (where  in  reality  less  than  M^  fringes 
would  be  crossed).  Figure  E-1  illustrates  the  geometry  of  the  calcu¬ 
lation  that  follows.  A  crossing  distance  variation  is  calculated 
first  and  then  a  crossing  time  variation  inferred. 


Figure  E-1  Measurement  Volume  Ellipsoid  Geometry. 
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Find  the  mean  vertical  crossing  distance ^  2*Y. 

_  a L  2  2  0  5  2 

Y  =  5  )  (L  -  (X/a)^  -  dZ*dX  /  (•n*a*L^/4) 

0  0 

=  ^  (Tr/4)*[L^  -  (X/a)^]*dX  /  (■n*a*L^/4) 

=  [L^*X  -  X^/(3*a^) I  ]  /  a*L^ 


=  2*L/3 


or  2*Y  =  4*L/3 


Find  the  mean  vertical  crossing  distance  squared,  ( 2  *Y ) 

-2  aL  /l^-  (X/a)^'  o  ,  ,  , 

Y  =  )  /  (L^  -  (X/a)^  -  Z^)  dZ*dX  /  (n*a*L^/4) 

0  0 

=  /  (2/3)*(L^  -  (X/a)^]  ^•^*dX  /  (Tr*a*L^/4) 

0 

9  Hi 

=  [8/(3*tT*a^*L^)  ]  *  (  J  (2/3)*[L^  -  (X/a)  ^  ^  *  ^*dX ) 

0 

=  [8/(3*TT*a^*L^)  ]  *  [3*Tr*a^*L'^/16] 


=  L^/2 


or  (2*y)^  =  2*L^ 


The  mean  crossing  distance  then  is  2*Y 
standard  deviation  of  crossing  distance  is: 


4*L/3  =  D.  The 


/2*L^  -  (4*L/3)^  =  0,471*L  =  0.354*D 


If  the  crossing  distance  has  a  standard  deviation  of  35%  about 
its  mean  value,  then  the  crossing  times  of  particles  will  have  a 
standard  deviation  about  their  mean  of  at  least  the  same  percentage. 
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The  distribution  about  T_(U)  is  an  error  in  the  weighting  function 

B  _ _ 

1/T„.  that  has  a  mean  of  0  and  a  standard  deviaition  of  0.35*T  (U). 

B 1  D 

The  error  is  independent  of  the  actual  velocity  U,  and  apparently 
increases  the  scatter  of  U  about  its  mean. 

tr;  .  [(Tu^  *  (.35*ir„)h°'^ 

=  1.06 
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APPENDIX  F:  FILTER  BIAS  ESTIMATES 


Assume  that  the  measurability  of  any 
the  signal  voltage  amplitude  after  filteri 
of  a  low  pass  filter  consider  figure  F-1 . 

Manufacturer  specifications  (TSI  I 
for  a  filter  roll  off  of  30  dB/octave  of 
every  factor  of  2  that  Fg^g  is  above  F^^ 
by  a  factor  of; 


velocity  is  proportional  to 
ng .  As  a  simplified  model 


nc.;  see  reference 
signal  voltage. 

, .  the  voltage  is 


3)  call 
So  for 
reduced 


R  =  10" dB/20)  ^  ^  Q32 


A  more  general  formula  is: 


R 


(F 


Sig 


/ 


F  .  \ 

Filter' 


-4.98 


(F-2) 


Measured  frequencies  ^^Sig^ 
component  magnitude  (U),  the  fringe 
frequency  shift  (FS)  : 


are  a 
spacing 


function  of  the  velocity 
(S) ,  and  the  component 


(F-3) 


Figure  F-1:  Low  Pace  Filler  Behavior  Idealization 


Substituting  equation  F-3  into  equation  F-2: 


(F-4a) 


(F-4b) 


(F-4c) 


Thus  for  a  given  random  sample  of  U  values,  the  best  estimate  of 
the  mean  is: 

N 

U  =II  U  /  N  (F-5) 

i  =  l 

while  the  arithmetic  mean  of  the  filter  biased  data  is: 

N  N 

b'  =11  R.ni.  /  Hr.  (f-6) 

i=l  i=l  ^ 


where  R. 


[  (U/S  +  FS)/F, 


-4 .98 


i  Iter 


FS  *  u^/s  > 

—  FS  +  ups  < 

- *-  otherwise 


Tht^  difference  between  the  two  mean  velocity  measurements  is  the 

filter  hi  IS  "error"  plotted  in  Figure  11.  The  curves  in  this  figure 

woro  rl( 1  oped  issuming  a  guassian  distribution  of  the  measured 

r  r«  a  !f  n('.7  aboi.it  its  mean  (mean;  F^ .  and  standard  deviation  CT, .  )  . 

Rig  Rig 


APPENDIX  H:  TRAVERSE  POSITIONING  ERRORS  DUE  TO  BENDING 

Positioning  error  values  P  ,  P  ,  and  P^  are  defined  as  the 

X  y  Z 

traverse  encoder  reading  minus  the  actual  LDV  measurement  location. 
These  errors  can  be  significant  and  in  a  worst  case  may  reach  25  mm 
or  more.  In  this  section  (as  previously),  the  positive  traverse 
coordinate  directions  are:  "X”  -  up;  "Y"  -  south,  on-axis;  "Z"  - 
east,  streamwise.  Careful  measurements  of  the  traverse  system 
structural  bending  were  made  with  a  very  accurate  inclinometer  (30 
seconds  accuracy) .  A  mill  table  and  calipers  were  also  used  to  meas¬ 
ure  directly  position  errors.  The  following  observations  were  made. 

1.  Over  the  range  of  possible  "X"  values  (-.25  m.  to  +.25 

m.)  there  is  no  evidence  of  structural  bending  changes  or 
therefore  direct  position  error  changes.  However  there  is  a 
clear  indication  of  an  indirect  effect  on  the  ^  2  “ 

P^,  error  relationships 

2.  Over  the  range  of  possible  "Y"  values  (-.50  m.  to  +.50 

m.)  the  Pj^  positioning  error  varies  from  0.0  to  +1.5  mm  and 
the  Py  error  varies  from  -7.0  to  +4.0  mm.  The  P^  position 
error  does  not  seem  to  change  significantly. 

3.  Over  the  range  of  possible  "Z"  values  (-.25  m.  to  +.25 

m.)  the  P^  positioning  error  varies  from  0.0  to  +1.5  mm  and 

the  P  error  varies  from  -12.0  to  +10.0  mm.  The  P.  position 
z  y 

error  does  not  seem  to  change  significantly. 

4.  Traverse  structure  inclination  changes  depend  on  where 
in  the  structure  you  measure  it.  The  structure  below  and 
including  the  inner  vertical  traverse  tube  bends  more  than 
the  outer  fixed  tube  and  structure  attaching  the  traverse 
system  to  the  carriage.  Inclination  variations  of  sections 
of  the  traverse  rails  are  slightly  larger  than  those 
measured  on  the  strut  or  probe  box. 

5.  There  is  indirect  evidence  of  P^ ,  P^ ,  and  P^  errors 
depending  on  the  pivot  position  or  distance  of  the  vertical 
traverse  tube  from  the  carriage  A-frame  mount.  In  pivot 
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position  3  (close  to  the  A-frame)  the  strut  inclination 

changed  by  12.5  minutes  over  the  entire  Y  traverse  range. 

In  pivot  position  11  (further  from  the  A-frame)  the  strut 
inclination  changed  by  15.0  minutes  over  the  entire  Y 
traverse  range.  Pivot  position  3  is  recommended  for  all  LDV 
tests  to  minimize  structural  bending  and  still  accommodate 
all  traverse  vertical  angles. 

6.  There  is  evidence  that  side  to  side  bending  angles  (9„e 

in  X-Y  plane)  are  substantially  independent  of  X  and  Z  posi¬ 
tion.  There  is  also  evidence  that  front  to  back  bending  an¬ 
gles  (0pg  in  X-Z  plane)  are  substantially  independent  of  X 
and  Y  position. 

7.  Both  positioning  error  and  bending  angle  measurements 

show  evidence  of  errors  due  to  hysterisis.  This  may  account 
for  much  of  the  data  scatter  and  represent  the  accuracy 
limit  of  the  following  correction  scheme. 

Simple  Model  of  Positioning  Error 

Assume  that  the  traverse  structure  is  completely  rigid  except 
for  bending  at  one  point,  a  "bending  point",  somewhere  along  the  axis 
of  the  vertical  traverse.  There  is  an  axis  or  Y,Z  position  where 
ther  are  no  moments  or  bending  about  this  point.  An  approximate 
calculation  and  symmetry  condition  place  this  "balanced  axis"  at  Y  = 
.1  m.  and  Z  =  0.0  m.  (see  figure  G-1).  The  sign  of  the  Y  location 
depends  on  whether  the  strut  is  north  (+)  or  south  (-)  of  the  meas¬ 
urement  point. 

Consider  for  each  traverse  location  that  a  move  is  made  from  the 
balanced  axis  to  that  traverse  location  with  no  bending  allowed. 
Then  the  structure  is  allowed  to  bend  (at  the  bending  point)  in 
response  to  tiie  moments  created  by  the  traverse.  All  movements  of 
the  LDV  measurement  volume  .aused  by  this  bending  are  the  positioning 
errors  at  that  traverse  location. 

Tn  this  model  bending  moments  and  hence  bending  angles  are  only 
functions  of  Y  and  Z  position.  X  position  contributes  only  to  the 


N 
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Positioning  Error  Predictions 

Since  measured  magnitudes  of  ^nd  0^^^  were  not  more  than 

0°  20',  the  following  approximations  are  very  good. 

(1  -  cos  =  (1  -  cos  9gg)  =  0.000  (G-2) 

sin  9p„  =  (in  radians) 

t  n  r  15 


sin  9^,, 


9 


(in  radians) 


*®FB 


(G-3a) 

(G-3b) 


(G-3c) 


p,  =  -  (Y  -  yo)*»ss  - 
Py  = 

=  -(X^  -  x)*e^3 
where  9gg  =  f(Y) 

Yq  =  0.201  m. 

OpB  =  g(z) 

Z  =  0 . 000  m. 
o 


There  are  three  missing  pieces  of  information  that  are  needed  to 
allow  this  model  to  predict  positioning  errors;  f(Y),  g(Z),  and  X^. 
These  can  be  empirically  deduced  from  a  set  of  data,  part  of  which 
has  already  appeared  in  figures  16  and  17. 

Based  on  optic's  platform  inclination  measurements  and  traverse 
encoder  readouts,  preliminary  functions  f'(Y)  and  g' (Z)  were  ob¬ 
tained.  A  least  square  error  fit  to  the  data  was  made,  assuming  a 
relational  form  of  angle  (0gg  or  9pg)  being  a  third  degree  polynomial 
of  position  (Y  or  Z  respectively) .  Figure  G-3  graphs  these  relation¬ 
ships  and  the  utilized  data.  These  optics  platform  inclination  an¬ 
gles  were  assumed  proportional  to  the  needed  inclination  angles  (9gg 
and  Qpg) •  The  following  multiplicative  constants  were  defined: 


^SS  “  ®SS  ^  ®SS 
^FB  ®FB  ^  ®FB 


The  unknowns  then  became  three  unknown  constants:  X^,  Kgg,  and  Kpg. 

The  main  data  set  included  positioning  error  measurements  as  a 

function  of  the  absolute  traverse  encoder  readouts.  It  was  obtained 

by  careful,  direct  measurements  utilizing  a  mill  table  and  calipers. 

P  and  P  resolutions  of  .1  mm  were  possible.  P„  resolutions  of 
y  z  X 

.25  mm  were  also  possible.  The  missing  constant  values  were  chosen 
as  those  that  minimized  the  sum  of  the  errors  squared  (difference 


is 


[i 


If 


ww 


WWg’JP 


between  measured  values  of 
values) . 


X  ' 


y' 


P^  and  the  formula  predicted 


Q  =  2.67  m. 

(G-4a) 

^  =  0.201  m. 

(G-4b) 

Q  =  O.COO  m. 

(G-4c) 

ss  =  1-2” 

(G-4d) 

SS  '  '^3*^  *  '‘2*’'  *  h**  *  '‘O 

(G-4e) 

,B  -  0.778 

(G-4f) 

FB  “  *  ®2*2^  * 

(G-4g) 

Implementing  Traverse  Positioning  Error  Corrections 

Accelerometers  have  been  permanently  attached  in  the  area  of  the 
probe  strut  to  optics  frame  connection.  These  can  be  used  to  monitor 
0gg  and  0pg  at  that  point.  The  initial  data  set  of  inclination  an¬ 
gles,  positions,  and  positioning  errors  has  already  been  displayed 
along  with  the  five  derived  constants,  and  two  functional  relation¬ 
ships  (eq.  G-4).  These  were  all  obtained  with  the  carriage  station¬ 
ary  and  the  strut  entirely  in  the  air. 

The  constant  values  (X 


,  Y  ,  Z„,  Kpjj,  and  Kl.„)  will  be  assumed 


O'  -o'  "o'  -'SS'  "IQ 

to  be  relatively  unchangeable  with  time  and  dynamic  strut  loadings. 

They  still  should  be  checked  every  couple  of  years  to  ensure  that 
traverse  wear  is  not  affecting  them.  This  is  also  true  of  the 
traverse  balance  point  derived  in  figure  G-1  (absolute  encoder 
readouts  Y  =  +/-0.1  m.  ;  Z  =  0.0  m.).  The  functions  for  0g^^  ,ind  0pg 
(figure  G-3)  will  only  serve  as  default  relationships  because  these 
functions  should  be  rederived  every  measurement  setup.  This  is  pos¬ 
sible  because  the  inputs  needed  to  derive  these  relationships  u. 
only  the  absolute  encoder  positions  and  the  -inqle  me  a  s' i  r  •■■th  ■  n  »  s  ‘  ■■■ 

permanently  attached  accelerometers.  These  inputs  ii  b.  r  -  i  •  : 

before  with  the  carriage  stationary.  Rut  the  bes’  -otr'  ‘I  'v.  . 

result  if  they  are  obta  ined  dur  ing  pre 1 im i n  a  t  y  ■  1 1  t  i  i  n  :  r  .  a . •  ■  *  • 

model  nhio  attached.  The  derived  functions  w i 1 .  '  .t-  ’  '  ■  ■  c 
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as  well  as  dynamic  effects. 

In  order  for  the  towing  carriage  software  to  do  the  positioning 
corrections  during  data  taking  the  parameters  of  this  section  must  be 
entered  into  the  computer.  The  "P”  selection  of  the  "U"  towing  tank 
menu  brings  up  either  a  "left"  or  "right"  parameter  screen  (Table 
G-1).  The  parameters  can  be  entered  directly  into  this  displayed 
screen  from  the  keyboard.  Note  that  each  entered  number  is  at  most  8 
digits  and  a  decimal  point  and  that  all  dimensions  are  in  inches. 
This  can  be  a  problem  with  with  the  3rd  Power  polynomial  coefficients 
which  can  be  on  the  order  of  10  .  So  in  Table  G-1,  all  polynomal 
coefficients  are  multiplied  by  1000  and  the  proportionality  constants 
are  divided  by  1000. 


G1  -  POSITION  CORRECTION  PARAMETER  SCREENS  (Right  =  North  Probe  Strut; 
Left  “  South  Probe  Strut) 


Traverse  Deflection  Parameters 
Side  =  Right 


Absolute  Position  of 
Relative  Home  in  (in) 
X  =  ?.???? 

Y  =  ?.???? 

Z  =  ?.???? 


Absolute  Position  of 
Balanced  Traverse  (in) 
X  =  105.2 
Y  =  7.91 
Z  =  0.0000 


Angle  Polynomial  Coefficients 
Y  -  Axis 

3rd  Power  =  .0001117 
2nd  Power  =  -.001377 
1st  Power  =  .04527 
0th  Power  =  -.1662 


Angle  Polynomial  Coefficients 
Z  -  Axis 

3rd  Power  =  .002172 
2nd  Power  =  .01154 
1st  Power  =  .4215 
0th  Power  =  0.0000 


Angle  Proportionality  Constants 
Y-Axis  =  .001273 
Z-Axis  =  -.000778 


Traverse  Deflection  Parameters 
Side  =  Left 


Absolute  Position  of 
Relative  Home  in  (in) 
X  =  ?.???? 

Y  =  ?.???? 

Z  =  ?.???? 


Absolute  Position  of 
Balanced  Traverse  (in) 
X  =  105.2 
Y  =  -7.91 
Z  =  0.0000 


Angle  Polynomial  Coefficients 
Y  -  Axis 

3rd  Power  =  .0001117 
2nd  Power  =  .001304 
1st  Power  =  .04468 
0th  Powe r  =  .16  50 


Angle  Polynomial  Coefficients 
Z  -  Axis 

3rd  Power  =  .002172 
2nd  Power  =  .01154 
1st  Power  =  .4215 
0th  Power  =  0.0000 


Angle  Proportionality  Constants 
Y-Axis  =  .001273 
Z-Axis  =  -.000778 
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EXPERIMENTAL  MEAN  VELOCITY  SCATTER 
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PROPELLER  INFLOW  SURVEY 


-  HIGH  SPEED  SHIP 

-  No  Propeller 

-  110  Data  Points 
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PROPELLER  INFLOW  SURVEY 

-  HIGH  SPEED  SHIP 

-  With  Propeller 

-  150  Data  Points 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


